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ABSTRACT 


The  performance  of  a  single-cylinder,  low  speed, 
spark  ignition,  internal  combustion  engine  has  been 
studied  using  lean  (i.e,  air-rich)  mixtures  of  propane 
as  the  fuelo  The  power  output  and  thermal  efficiencies 
have  been  determined  at  various  compression  ratios  and 
fuel-air  ratios .  A  comparison  is  also  made  with  the  en¬ 
gine  when  burning  gasoline. 

Maximum  operating  fuel  economy  is  obtained  at  a 
fuel-air  ratio  of  0o04  lb.  of  propane  per  lb.  of  air, 
regardless  of  the  compression  ratio.  It  is  also  shown 
that  the  overall  performance  of  an  engine  may  be  im¬ 
proved  by  burning  propane  as  the  fuel  at  a  higher  com¬ 
pression  ratio  than  burning  gasoline  at  a  lower  CR. 

A  theoretical  analysis  is  also  shown  for  obtain¬ 
ing  "cycle"  temperatures,  indicated  thermal  efficiency 
and  brake  mean  effective  pressure. 
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CHAPTER  I 


INTRODUCTION 

The  fundamental  problem  considered  in  this  investi¬ 
gation  is  the  experimental  determination  of  the  effects 
of  lean  fuel-air  ratios  of  gaseous  propane  (C«,H_)  when 
used  as  a  fuel  in  a  spark  ignition  four-cycle  internal 
combustion  engine 0 

lol  HISTORY 

The  spark  ignition  four-cycle  reciprocating 
internal  combustion  engine  cycle  is  usually  referred  to 
as  the  Otto  cycle,  after  N0A0  Otto  who  is  believed  to 
have  made  the  first  successful  internal  combustion  engine 
in  about  1876 0 

A  detailed  explanation  of  the  ideal  Otto  cycle 
and  deviations  from  the  ideal  cycle  are  given  in  Chapter  II , 


lol-X  Theoretical  Cycles 
Goodenough  and  Baker 


improved  the  understand¬ 


ing  of  internal-combustion  engine  analysis  by  considering 
actual  mixtures  of  gasoline,  clearance  gases,  variable 


*  Numbers  in  square  parenthesis  designate  References 


on  page  56  . 


. 
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specific  heats  and  chemical  equilibrium.  Keenan  and  Kaye 
prepared  a  set  of  highly  accurate  gas  tables  giving 
the  thermodynamic  properties  of  air  and  the  combustion 
products  of  octane  for  fuel-air  ratios  of  25  and  50  percent 
of  air  for  complete  combustion  at  low  temperatures  i,e, 
they  did  not  consider  the  dissociation  of  CO^  and  H^O, 


Hottel,  Williams  and  Satterfield 


presented  a  series 


of  charts  for  hydrocarbon -air  combustion.  Charts  are 
given  for  fuel  expressed  as  for  six  fuel-air  ratios 

ranging  from  80  percent  to  150  percent  required  fuel, 

A  theoretical  analysis ,  using  propane  as  the  fuel 
is  presented  in  Chapter  II. 


1,1-2  Experimental  Work. 

During  the  past  several  decades  numerous  experimental 
analyses  have  been  conducted  on  spark-ignition  internal 
combustion  engines,  using  various  fuels,  but  mainly  gaso¬ 
line,  the  most  commonly  used  fuel  in  a  practical  engine. 

The  prime  objective  of  any  engine  is  power  output  and 
since  this  is  in  direct  proportion  to  the  number  of  cylind¬ 
ers,  mean  effective  pressures,  length  of  stroke,  area  of 
cylinder  and  number  of  revolutions  per  minute,  an  increase 
in  any  one  of  these  for  a  particular  fuel  will  result  in 
a  power  increase.  Practical  considerations  however  limit 
the  number  of  cylinders  and  size  of  cylinders.  Increasing 


3 


the  speed  will  increase  the  power  output,  up  to  an  op¬ 
timum  value.  However,  prior  to  this  the  effect  of  speed 
increase  has  been  to  cause  the  efficiency  to  be  reduced,, 
The  greatest  possibility  for  increase  of  horsepower  per 
unit  of  engine  weight  lies  in  increasing  the  mean  effect¬ 
ive  pressure. 

The  most  obvious  way  to  increase  the  mean  effective 
pressure  is  to  increase  the  compression  ratio  of  the 
engine,.  For  a  particular  fuel,  however,  there  is  a  max¬ 
imum  compression  ratio  above  which  the  fuel  will  detonate 
(Engine  Knock) 0  Detonation  results  in  a  loss  of  power 
and  if  severe,  the  engine  will  be  damaged,.  Many  factors 
determine  when  a  particular  fuel  will  knock,,  Some  of 
these  factors  are  combustion  chamber  design,  engine  speed, 
spark  timing,  cylinder  temperature,  air-fuel  ratio  and 
air-fuel  mixture  temperature „  Gasoline  with  knock  in¬ 
hibitors  are  now  available  which  will  allow  engines  to 
operate  without  knock  at  a  higher  compression  ratio  than 
if  the  gasoline  were  used  without  the  inhibitor „  Many 
investigations  have  been  conducted  regarding  Engine 


Knock 


4,  5,  6,  7 


Propane  has  a  higher  detonation  resistance  than 
all  but  the  most  expensive  aviation  gasoline.  Propane 
also  has  the  other  characteristics  necessary  for  a  fuel 
in  an  I.C.  engine  i.e.  adequate  volatility  without  pre¬ 
heating,  low  sulphur  content,  has  high  heating  value. 


...  •:  •  '  *C.  '  '  R  • 
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readily  available  at  an  economical  cost,  and  the  added 
advantage  over  gasoline  in  that  its  combustion  is  more 


complete 


8 


Thus  carbon  deposits  in  the  combustion 
chamber,  and  on  the  valves  and  piston  head  are  notably 
reduced  resulting  in  less  dirty  engine  oil,  resulting 
in  fewer  oil  changes  and  less  engine  wear0  There  is  also 
a  reduction  in  the  noxious  exhaust  odors,  smoke  and  soot0 
Another  advantage  of  propane  is  that  it  has  a  slow  rate 
of  flame  travel  during  combustion,  resulting  in  a  de¬ 
crease  in  peak  pressure  but  resulting  in  an  increase  in 
mean  effective  pressure,. 

Operating  variables  such  as  spark  advance,  exhaust 
pressure,  inlet  pressure,  inlet  mixture  temperature  also 
have  an  influence  on  the  mean  effective  pressure  and  the 
engine  efficiency 0  A  very  thorough  comparison  of  the  ef¬ 
fects  of  these  variables  is  shown  in  Taylor 

The  heating  value  of  a  fuel  is  not  in  direct  pro¬ 
portion  to  the  power  output  which  can  be  obtained  from  the 
use  of  that  fuel  in  a  particular  engine  because  of  other 
factors  which  affect  the  power  output.  The  heating  value 
is,  however,  a  direct  measure  of  the  quantity  of  fuel 
required  the  higher  the  heat  value  the  less  fuel  needed 
to  do  the  same  work.  Propane  has  a  slightly  higher  heat¬ 
ing  value  per  pound  than  gasoline  but  a  lower  value  if 
compared  on  a  volume  basis. 


. 


■ 


£ 


Adams  and  Boldt 


10 


conducted  extensive  tests  on 


three  industrial  engines /  including  one  engine  at  three 
compression  ratios  using  gasoline  and  liquified  petroleum 
gas  (Butane-propane) ,  They  found  that  at  a  given  com¬ 
pression,  ratio ,  the  engines  developed  3-1/2  -  5%  less 
power  at  high  speeds  on  L0Pc-gas  than  on  gasoline.  How¬ 
ever  compared  to  gasoline  the  L„P,-gas  reduced  brake 
specific  fuel  consumption  on  a  mass  basis  up  to  12%  at 
low  speeds.  At  high  speeds  fuel  consumption  was  reduced 
by  0  -  9%,  On  a  volume  basis,  the  L.P,  gas  consumption 
increased  due  to  the  lower  specific  weight  of  the  L„P,-gas, 
They  also  found  that  increasing  the  compression  ratio  from 
7,5  to  11,5  increased  the  power  12%  and  reduced  the  brake 
specific  fuel  consumption  by  11%. 

Because  of  the  certain  advantages  of  propane  as  a 
fuel  for  internal  combustion  engines  it  is  the  purpose  of 
this  thesis  to  obtain  data  on  an  engine  operating  at 
maximum  economy  fuel-air  ratio.  Although  this  fuel-air 
ratio  is  at  a  lower  than  maximum  power  such  information 
is  necessary  and  of  importance  for  vehicles  when  operating 
at  cruising  speed  and  is  also  very  important  to  the  designer 
and  owner  of  stationary  engines  wishing  to  operate  at 
maximum  economy  over  an  extended  period  of  time. 


emulov 

oni  :ii; 


CHAPTER  11 


THEORETICAL  ANALYSIS 


2.1  GOVERNING  EQUATIONS 
2.1-1  Otto  Air-Standard  Cycle 

The  simplest  and  crudest  analysis  of  the  Otto  cycle 
considers  the  working  medium  as  air  or  at  least  a  gas 
having  the  properties  of  air.  It  is  further  assumed  that 
the  air  has  a  constant  specific  heat. 

This  cycle  is  shown  on  the  P-V  and  T-S  diagrams  of 
Fig.  2.1. 


T 


PIG.  2.1  OTTO  CYCLE 

The  medium  is  compressed  reversibly  and  adiabatically ,  i.e, 
isentropically ,  process  1-2.  Heat  is  added  at  constant 
volume,  process  2-3,  while  the  piston  is  momentarily  at 
rest  at  "top-dead  center".  Ignition  occurs  prior  to 
point  3  in  the  actual  engine.  Isentropic  expansion  takes 
place,  process  3-4,  and  process  4-1  is  the  rejection  of 
heat  while  the  Pi  ston  is  returning  to  its  initial  position 
at  "bottom-dead  center". 
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From  the  First  Law  of  Thermodynamics  for  a  closed 
system  and  assuming  no  changes  in  kinetic,  potential  or 
chemical  energies  and  assuming  adiabatic  compression  and 
expansion  the  net  work  of  the  system  may  be  obtained  from: 


Net  W  =  Q.  -  Q  . 

in  out 


where 


in 


me  (T0  -  T0)  and  Q  ,  =  me  (T.  -  T,  ) 

v  3  2  out  v  4  1 


The  thermal  efficiency  ( t]) 


net 


W 


therefore 


n  = 


°in  Qout 

Qin 


out. 

Qin 


(2.1) 


=  1  - 


mCv(T4 


Tl> 


me  (T 
v  3 


T  ) 
1  2 


so 


=  1  - 


T  -  T 

4  1 

T  -  T 

1 3  2 


(2.2) 


x.h 

T 

2 


^/Ti 


-  i 


-  1 


But  from  isentropic  relationships 
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so 


1 


T 


T 


2 


1 


1 


V 


1 


But 


V 


2 


rr~  -  CR  (Compression  ratio)  . 
V3 


1  -  CR 


1-k 


(2.3) 


Thus  in  the  air- standard  Otto  cycle  the  thermal 
efficiency  is  a  function  of  only  the  compression  ratio  and 
the  efficiency  may  be  increased  by  increasing  the  compression 
ratio.  Although  the  thermal  efficiency  will  increase  to  a 
maximum  value  of  one  as  CR  increases, the  rate  of  increase 
of  T)  decreases  as  CR  increases,  and  any  increase  is  prac¬ 
tically  insignificant  as  the  CR  approaches  15  or  20. 

2.1-2  Otto  Fuel-Air  Cycle 

The  spark-ignition  engine  deviates  from  the  air-stand¬ 
ard  cycle  in  many  ways.  Some  of  the  more  important  ways 
are  as  follows? 

1.  The  medium  consists  of  a  mixture  of  fuel,  air  and 
residual  gases. 

2.  The  specific  heats  of  the  individual  gases  are  dif¬ 
ferent  and  vary  with  temperature. 

3.  The  combustion  process  replaces  the  heat  transfer 
process.  Combustion  takes  time  and  is  thus  initiated 


before  top-dead  center. 


- 
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4.  After  combustion ,  i„e.  during  adiabatic  expansion 
and  exhaust,  the  medium  is  no  longer  a  mixture  of 
fuel  and  air  but  air  and  products  of  combustion » 

50  The  reaction  may  not  be  complete „ 

The  fuel-air  cycle  is  an  idealized  thermodynamic  pro¬ 
cess  resembling  what  is  occurring  in  an  ''ideal'1  engine  be¬ 
cause  in  a  real  engine  the  process  is  not  cyclic  in  the 
thermodynamic  sense „ 

Further  deviations  take  place  in  the  actual  engine 
due  to s 

10  The  inlet  and  exhaust  processes  will  result  in  a 
pressure  drop  through  the  valves  thus  resulting  in 
a  loss  of  worko 

2,  The  piston  must  do  work  on  the  mixture  to  get  it 
out  of  the  cylinder  and  this  is  more  than  the  work 
done  in  the  cylinder  by  the  mixture  at  intake,  re¬ 
sulting  in  a  net  loss  of  worko 

30  Heat  transfer  is  involved,  so  the  compression  and 
expansion  processes  are  not  isentropic,  as  assumed. 

Fig,  2.2  illustrates  the  deviation  of  the  actual  cycle 
from  the  ideal  cycle  as  it  affects  the  P-V  diagram. 

2.2  ANALYSIS  USING  PROPANE 


The  reaction  equation  for  propane  reacting  with  the 


r 


10 


V 


IG 


2  o  2 


ACTUAL  FUEL-AIR  CYCLE  LOSSES 


4 


11 


correct  amount  of  air  for  complete  combustion  i.e. 
stoichiometric  air-fuel  ratio,  iss 

C3H8  +  5  02  +  3o76  (5)N2  3C02  +  +  3.76(5)N2  (2  „  2-1) 

(assuming  air  21%  oxygen  and  79%  nitrogen  by  volume,  where 
the  %  nitrogen  includes  the  inert  gases  in  the  air) .  If 
there  is  an  excess  of  air,  oxygen  will  appear  on  the  right 
hand  side  of  the  equation  and  the  nitrogen  on  both  sides 
will  also  be  increased. 

If  there  is  insufficient  air  for  complete  combustion 

i„e.  a  rich  fuel-air  mixture  the  number  of  moles  of  CO^ 

will  be  reduced  and  there  will  be  CO  in  the  products  and 

a  reduction  in  the  number  of  moles  of  N_ , 

2 

In  an  engine,  even  though  operating  on  an  ideal  cycle, 
there  will  be  some  gases  left  in  the  cylinder  from  the 
previous  cycle.  These  gases  CO^ ,  H^O ,  N2  and  CO  or  Q^f 
depending  on  whether  a  rich  or  lean  mixture  is  being 
burned,  are  at  a  higher  temperature  (T^)  and  although  the 
volume  is  small  (approximately  3%  of  V^)  there  is  an  ap¬ 
preciable  increase  in  energy  and  temperature  of  the  in¬ 
coming  mixture  when  mixed  with  the  residual  gases. 

Assuming  100%  air  and  "f"  as  the  fraction  of  residual 
gases  remaining  in  the  cylinder  the  reaction  equation  2.2-1 


now  becomes: 
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C3Hq  +  5  02  +  18,8N2  +  £(3CQ2  +  4H20  +  18. 8N2) 

3(1  +  f)C02  +  4(1  +  £)H20  +  18.8(1  +  £)N?  (2.2-2) 


The  temperature  of  the  incoming  fuel -air  mixture  may  be 
measured  and  thus  the  enthalpy  obtained.  However,,  the 
temperature  of  the  residual  gases  and  "f"  have  to  be  esti¬ 
mated  and  later  confirmed,  i.e.  a  trial  solution. 

The  internal  energy  at  T=^  may  be  expressed  ass 


H  ,  +  U_ 

charge  Res.  gas 


MRT 

but  PV  =  MRT  .  °  .  P_  =  -7“ 

I 


and  V  “ 

o  2 


o 


®  o 


M (1545) 

778 


M (1 o 986)T1 


=  M  (1  o  986)T 


where  CR  is  the  compression  ratio. 


o  e 


U 


1 


“fuel  +  H02  +  HN2  +  “res . gas  “  **(!. 986)^1  -  - 


(2.2-3 


At  an  assumed  temperature  T^,  can  thus  be  obtained.  Then 


considering  the  entire  mixture  at  T^,the  sum  of  the  internal 


energies  of  the  gases  can  be  obtained  froms 


U 


1 


tL 


u  TT  ^  IT  TT 

fuel  02  N2  C02  "  ^H20 


This  process  is  repeated  at  various  T, "s  until  both  cal- 

l 


culated  1 s  are  of  the  same  value.  The  value  of  T  thus 


obtained  is  taken  as  T. 


Xsentropic  compression  from  T^  to  T2  gives  a  means  of 


obtaining  T 2  <, 


T, 


V 


1 


0 


MC  In  ■=- 
v  T. 


MR  In 


2 


V, 


T, 


V, 


V. 


1.8 


M  c  In  -r 
v  T, 


MR  In 


V. 


MR  In 


V, 


MR  In  CR 


Various  values  of  T2  are  assumed  until 


T 


Me  In  _ 
v  T 


2 


1 


M  R  In  CR 


2o2-< 


The  energy  at  T'2  is  then  calculated  from 


U 


u  +u  +u  +  u  +u 

c3h8  02  N2  C02  h20 
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U  and  S  are  calculated  using  Cp  and  Cv  as  a  function  of 


temperature 


11 


or  from  the  tables  based  on  Cp  and  Cv 


as  functions  of  temperature 


12 


,  (Table  datum  is  520  R) 


The  process  2-3  of  Fig.  2*1  represents  the  combustion 
of  the  fuel  at  constant  volume*  In  the  ideal  case  this 
process  is  assumed  to  be  adiabatic  and  the  gas  mixture 
is  in  equilibrium  when  state  3  is  reached* 

The  mixture  at  the  high  temperature  of  T^  will  con¬ 
tain  unburned  CO  and  H0  i.e.  the  CO^  and  H^O  dissociates 
into  CO  and  and  and  ,  respectively. 

If  x  and  y  represent  the  degrees  of  dissociation  of 
CO^  and  H^O,  respectively,  the  reaction  equation,  2.2-2 
can  be  written  ass 


C~H0  +  50o  +  18. 8N«  +  f ( 300 «  +  4H„Q  +  18.8FL) 
3  o  2  2  2  2  2 


3 (1+f) (1-x)  C02  +  3 (1+f ) (x) CO  +  4(l+f) (l-y)H20  +  4(l+f)(y)H2 

+  (1+f)  ^  Q2  +  (1+f)  (18o8)N2  (2.2-5) 

The  criteria  for  chemical  equilibrium  is  that  the 
Gibbs  function  (Gibbs  free  energy)  for  the  system  must  be 
a  minimum.  The  Gibbs  function  is  defined  ass 

G  ~  H  -  TS 

By  differentiating  and  substituting  it  can  be  shown  that 
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dG  =  VdP  -  SdT 


At  the  final  approach  to  equilibrium  the  process  may  be  con¬ 
sidered  isothermal,  i „ e „  dT  =  0. 

o 

dG  =  VdP 

RT 

Assuming  a  perfect  gas  relationships  V  =  -p- 

8  SP 

dG  —  RT 

P 

Integrating  this  equation  gives? 

G2  “  G1  “  AG  —  RT  In  (P2/P  )  ° 

For  the  reaction  aA  +  bB  cG  +  dD 


AG 


2 


AG^  -  AG^  +  AGq  -  -  AGg 

=  CRT  In  (P2/P1)c  +  dRT  In  (P2/Pi)d 
-  aRT  In  (P2/P1)A  -  bRT  In 


RT  In 


'P2C^  ^P2D^ 

‘P2A^  ^P2B^ 


RT  In 


>  ) c  (p  ) d 

lCj  v  ID' 


P1A^ 


IB' 


Arbitrary  values  may  be  assigned  so  that  the 
last  term  has  a  value  of  zero. 


i.e.  Plc  =  P1D  -  P1A  =  P1B  =  1  atmosphere. 
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AG2  -  AGX 


RT  In 


(P2C)  ^P2D^ 

(p2A)  (P2E^ 


Considering  point  2  as  the  condition  at  equilibrium 

0 

where  there  is  no  change  in  free  energy,  AG2  becomes  zero 
and  AG^  becomes 


-AG 


0 


RT  In 


(pA)a  (PB)b 


or 


0 


RT 


(pc)°  (pD)d 


V 


b 


B' 


R 


P 


where  represents  the  pressure  equilibrium  constant » 
Since  the  partial  pressures  of  a  mixture  of  gases  may  be 
assumed  proportional  to  the  mole  fractions 


K 


d 


(A)  (B) 


b 


;2„2-6) 


K  is  a  function  of  temperature  and  values  of  K  for  various 
P  —  P 


gases  are  tabulated 


13 


Since  the  deviation  of  the  actual 
gas  from  the  ideal  gas  is  slight,  the  compressibility  cor¬ 
rection  factor  Z  may  be  taken  as  1  in  the  range  of  tempera¬ 
tures  and  pressures  encountered  in  combustion. 


K 


K 


..  '  I  1  . 
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0 


2' 


Considering  the  dissociation  of  200^  into  2C0  and 


CO 

2 

°2 

CO 

2 

(2.2-7) 


Substituting  for  the  values  of  the  number  of  moles ,  from 
equation  2,2-5,  equation  2,2-7  becomes: 


K 


CO, 


3  (1+f)  (x)P3l 

2 

"  3 (1+f) x  +  4  (1+f) y 

M3  J 

L  2m3  J 

3 (1+f) (1-x) 


M. 


x 


1-x 


|  (l+£)x  +  |  (1+f)  y 


M, 


but 


P. 

PL 


RT  P  V 

where  R  «  and 

V3  M2T2  2  3 


O  © 


K 


CO, 


1 2  r 


1-X 


~  (l+f)x  +  “  (l+f)y 


P2T3 

M2T2 


(2,2-8) 


The  equilibrium  constant  for  the  dissociation  of  2HL>0  may 
be  expressed  in  the  same  manner. 

For  reasons  to  be  explained  later  it  is  advantageous 
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to  determine  the  water-gas  equilibrium  constant  K^_,  rather 
than  q.  Considering  the  two  equations % 


2  CO.  2C0  +  0 

2  2 

2H20  2H2  +  02 

Subtraction  gives 

2CQ2  -  2H20  2C0  -  2H2 

or  H.  +  C0o  HO  +  CO 

2  2  2 


In  terms  of  x  and  y  and  mole  fractions  this  equation  becomes 


a  iLlizXL 
(l-x)y 


(2.2-10) 


Values  of  at  various  temperatures  are  also  tabu¬ 


lated 


13 


Theoretically ,  the  „  Q 2  and  N2  will  dissociate  into 
their  respective  atoms  and  small  amounts  of  OH,  NO,  NH^  and 
CH^  may  form,  but  at  high  pressures  as  found  in  the  Otto 
cycle  f  this  dissociation  is  very  slight  and  is  neglected 
in  this  analysis. 

may  be  obtained  by  using  the  first  Law  of  Thermodynamics, 
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expressed  as 


U,.  +  CE-  , 
2  fuel 


U. 


and  the  CE  of  the  fuel  are  known.  The  CE  is  ob¬ 


tained  from  CE  +  TJ  +  U 


a 


fuel 


U  +  U  +  Heating 

\**\J  ^  iri  p 


Value  of  fuel . 


Heating  Value  of 


878,749  BTU/mole 
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U. 


UC02  +  uco  +  uh20  +  UH2  +  u02  +  \ 


3(l+f)  (1-x)  uco  +  3 


x  uco  + 


4  (1+f ) y  uH  +  (1+f  ■  +  4^)  uQ  +  18.8 


(1+f ) u 


N, 


(2.2-11; 


T^  ,  x  and  y  are  yet  unknown,  therefore,  this  equation 
cannot  be  solved. 

Considering  equations  2  <,2-8,  2  <,2-10  and  2.2-11  there 
are  three  equations  with  five  unknowns  ,  x,  y  and 
and  Kcq  ,  since  the  8,u’s"  of  equation  2.2-11  are  all  known 


2 


functions  of  temperature.  A  trial  solution  is  involved 
to  obtain  .  T^  is  assumed,  and  the  internal  energies  of 
the  various  gases  in  equation  2.2-11  are  obtained  by  using 
the  equations  u  - 


erature 


11,12 


=  C  dT  where  C  is  a  function  of  temp- 
J  P  P 

The  x,  y  and  numerical  terms  are  added 


algebraically  resulting  in  an  equation  of  the  form 
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°2  +  CEC3H8 


^terras  x  /  xterms  +  ^  /  ^terms 


Therefore  y  may  be  obtained  in  terms  of  x 

Number  -  x  \  x. 


i  o  e , 


y 


terms 


2-12) 


terms 


From  equation  2,2-10,  y  can  also  be  obtained  in  terms  of  x 
as  can  now  be  obtained  for  the  assumed  value  of  . 


x 


^G  X^KWG 


!2.2-13) 


Equating  equations  2,2-12  and  2,2-13  eliminates  y  and 
results  in  a  quadratic  equation  in  terms  of  x  from  which 
a  value  of  x  may  be  obtained,  y  may  then  be  calculated 
from  either  equation  2,2-12  or  equation  2,2-13,  Using 
these  values  of  x  and  y  in  equation  2,2-8  a  value  is  ob¬ 


tained  for  K 


CO, 


This  procedure  is  repeated  at  various 


H‘  0  s  until  K 


as  tabulated 


C?2 

13 


calculated  is  approximately  equal  to  K 


CO, 


When  various  ' s  are  obtained,  the 


correct  T^ ,  x  and  y  may  be  obtained  by  linear  interpola¬ 
tion  or  by  graphical  means ,  such  that  calculated 

equals  Kr  from  the  tables, 

CO  2 

Process  3-4  is  assumed  to  be  isentropic  expansion, 
therefore  -  S^,  Knowing  T^ ,  and  y^  the  entropy  at  3 
may  be  calculated  by  the  equations 


tr nr  ♦  3  'cT  tm  ' 


. 
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s. 


0  _ 


MR  In  pp„  + 


McoCEco 


T 


+ 


3 


M  CE 

tl ^  a 


:202-14) 


where  0  is  M 


C  dt 
_J£ _ 

T 


and  pp  is  the  partial  pressure  of 


the  constituent  gas. 

The  determination  of  again  involves  a  trial  solution 
since  the  degrees  of  dissociation  of  the  C0«  and  Ho0  cannot 

Z  Z 

be  determined  until  after  is  obtained*  Assuming  a  value 

of  T^,  Kco  and  K^G  may  be  obtained  from  the  tables*  Using 
2 

equation  2*2-13,  y  may  be  obtained  in  terms  of  x*  Substi¬ 
tuting  in  equation  2*2-8  results  in 


x 


2 


CO 


2 


1 —x 


3  2  (l+f)x 

.2 (1+f)x  *  - 15 x. 


p 

M 


4 


but 


P, 

M~ 


RT 


V 


4 


RTZ 


¥1 

MiTx 


T 

V, 


4 


t4(i 


M1T1 


since  P.  “1*  atm. 
JL 


therefore 


K 


CO. 


'  \2  r-  y  (l+f)x 

T~)  f  (i+f)x  +  - - ir-TZ - r 

.  x/  Lz  ^WG  '  wG 


lx 


T 


MiTi 


(2.2-15) 


Various  values  of  x  are  assumed  and  substituted  until 


as  calculated  equals  the  value  of  K  obtained  from 


CO. 


CO 


2 
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the  tables 0 

Using  the  assumed  value  of  and  y4  as  obtained,, 

the  entropy  at  may  be  calculated  by  the  equation 


<t>  _ 


-r  In  r  r  ■  Mco  CEC0  ,  MfI2  CEh2 
MR  In  pp  4*  - — -  + 


T 


T. 


If  exceeds  S^,  is  assumed  lower,  and  vice  versa. 

The  above  procedure  is  repeated  until  S^  “  S^,  thus  giv¬ 
ing  the  temperature  T^, 

At  4,  the  end  of  the  expansion  process,  the  exhaust 
valve  opens  and  the  combustion  products  escape  from  the 
cylinder  and  the  pressure  decreases  to  atmospheric  pres¬ 
sure,  The  gases  in  the  cylinder  expand  isentropically 
to  atmospheric  pressure.  Some  of  these  gases  remain  in 
the  clearance  space  in  the  cylinder.  Considering  this  as 
condition  5  it  is  possible  to  calculate  V^, and  knowing 
check  on  the  assumed  value  of  "f”  and  also  the  assumed 

value  of  Tc.  Since  Tr  is  in  the  order  of  2200 °R  it  is 
5  5 

reasonable  to  assume  the  amount  of  dissociation  of  CC>2 
and  H^O  is  negligible. 

For  an  assumed  T,_  the  entropy  of  5  may  be  calculated 

using 


S 


5 


MR  In  pp. 


Various  T  '  s  are  assumed  until  a  value  is  obtained  which 
o 


- 


results  in  S_  -  S~0 

5  3 

The  fraction  of  residual  gas,  f 
obtained 


V5  ' 


V 


2 


CR 


M.RT, 
1  1 


140  7  (144) CR 


may  now  be 


and 


V, 


MCRTC 

D  D 


14.7 (144) 


therefore 


V, 

Vr 


2 


W 


CR  (T 


5  5 


2-16' 


If  Tj.  and  "fdiffer  greatly  from  the  initially  assumed 

values  of  "flsand  Tr  the  entire  calculations  must  be  re- 

5 

peated  until  there  are  negligible  changes  in  T'and  T^o 
The  thermal  efficiency  is  defined  as 


__  Net  Work 
^  Heat  Added 


3 


U4> 


r2  -  ur 


Heat  Value  of  Fuel 


(2,2-17) 


The  heat  value  of  propane  is  taken  as  the  lower  heat 


value  of  19929 


14 


Knowing  the  temperatures  T^ , 


T^ ,  and  the  internal  energies  may  be  calculated  and 
the  thermal  efficiencies  obtained  by  equation  2 „ 2-17, 

The  mean  effective  pressure  may  be  obtained  from 


net  work 
displacement 


m  o  e  o  p  o 
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(U3  "  V  "  (U2  '  V 

therefore  m.e.p.  =  - - - 

(V1  "  V2}  144 


(2 . 2-16) 


A  sample  calculation  is  given  in  Appendix  A. 

A  computer  programme  in  Fortran  IV  language  was  written 
to  carry  out  the  above  calculations  for  various  percentages 
of  air,  compression  ratios,  and  various  fuels 0  The  computa¬ 
tions  were  done  on  the  IBM  7040-1401  system  in  the  Depart¬ 
ment  of  Computing  Science,  University  of  Alberta*  Appendix 
B  shows  this  program*  The  results  are  plotted  in  Chapter 
IV*  Iterations  were  continued  until  there  was  less  than 

O 

a  1  temperature  variation  in  any  particular  assumed  temp¬ 
erature  or  ^f/o  variation  in  any  other  assumed  value. 

Convergence  is  illustrated  by  means  of  the  table  which 
follows.  These  results  are  at  a  compression  ratio  of  10. 
Results  at  other  compression  ratios  showed  similar  effects. 


Percent  air 

Thermal  Efficiency 

20%  air  Intervals  10%  air  Intervals 

1%  variation 

^2%  variation 

100 

36,8 

36.9 

36.9 

110 

3  9  *  5 

39.5 

120 

43.7 

41.8 

41.8 

130 

43.6 

43.6 

140 

58.3 

44.8 

44.8 

150 

45.0 

45.0 

160 

46.6 

46 . 6 

46.6 

170 

46.2 

46.2 

180 

46.5 

46.9 

46.9 

190 

47.7 

47.7 

200 

46.2 

47.9 

47,9 

' 


CHAPTER  III 


EXPERIMENTAL  WORK 


3  o  1  EQUIPMENT 

The  engine  used  in  conducting  these  tests  was  an 
""ASTM-CFR-48  Engine"  located  in  the  Mechanical  Engineering 
Laboratory  of  the  University  of  Alberta.  The  engine  is 
connected  to  a  series  wound  D.C.  motor  capable  of  starting 
the  engine  and  absorbing  the  power  developed  by  it  (dyna¬ 
mometer)  .  The  engine  has  a  bore  of  3.25  in. ,  stroke  of 
4.50  in.,  variable  spark  adjustment  and  a  variable,  from 

4  to  10,  compression  ratio.  The  standard  CFR  engine  cool¬ 
ing  system,  evaporative  cooling  was  used  to  maintain  the 
coolant  temperature  at  208°F  (See  Fig.  3-1) . 

The  air  flow  into  the  air  inlet  surge  tank  was 
measured  by  means  of  an  orifice  and  inclined  manometer. 

Gaseous  propane  was  introduced  into  the  air  stream, 
at  the  throat  of  the  carburetor f  to  ensure  adequate  mixing. 
The  propane  flow  was  measured  by  means  of  a  "Precision" 

Wet  Test  Gas  Meter,  manufactured  by  Precision  Scientific 
Co. ,  placed  in  the  propane  gas  line  after  the  pressure 
reducer  and  before  the  engine  carburetor.  The  meter  cap¬ 
acity  was  0.250  cu.  ft.  per  revolution  of  the  sweep  dial, 
with  a  maximum  pressure  of  8  in.  of  water. 

The  combustion  chamber  mean  temperature  was  obtained 


■ 
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by  a  chrome  1  -a lumel  thermocouple,  placed  in  the  top  of 
the  cylinder,  at  the  "knock -meter  connection".  The  thermo¬ 
couple  readings  were  obtained  on  a  Leeds  and  Northrup 
Millivolt  potentiometer o 

The  time  of  runs,  for  the  combustion  of  0,500  cu, 
ft,  of  gas,  and  the  number  of  revolutions  per  run  were 
obtained  on  the  timing  unit  on  the  control  panel  (see 
Fig,  3-2)  , 

Exhaust  gas  analyses  were  obtained  by  means  of  a 
Hays  Gas  Analyzer  (Orsat)  ,  capable  of  determining  the 
percent  by  volume,  of  CO^ ,  O ^  and  CO, 

3,2  EXPERIMENTAL  METHODS 

Initial  tests  were  conducted  with  the  engine  as  pro¬ 
vided,  Later  tests  were  done  with  the  engine  exhaust 
valve  replaced  by  a  mild  steel  exhaust  valve,  manufactured 
locally,  to  determine  the  effects  of  oxidation,  if  any. 

No  apparent  changes  occurred  in  the  operation  of  the  engine 
as  a  result  of  the  changing  of  the  exhaust  valve,  A  photo 
of  this  valve  showing  the  thermocouple  wire  is  shown  in 
Fig,  3,3, 

3,2-1  Test  Procedure 

After  preliminary  testing,  the  engine  oil  was  re¬ 
newed  and  the  micrometer  gage,  indicating  the  compression 
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Front  View 


■Side  View 


FIG,  3-1 


CFR -ENGINE  AND  TEST  EQUIPMENT 
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PIG  c  3-2  ‘CONTROL  PANEL 
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ratio ,  was  calibrated 

Prior  to  a  test  run  the  engine  was  operated  for  a 
considerable  time  to  ensure  it  was  at  operating  tempera¬ 
ture,  The  compression  ratio  was  set  at  the  desired 
value  and  the  gas  pressure  regulator  adjusted  to  give 
the  desired  rate  of  fuel  flow.  The  load  was  then  adjusted 
to  give  the  desired  r,p,m.  (600  for  all  tests).  The 
spark  was  adjusted  to  give  the  maximum  power  output,  indi¬ 
cated  by  the  lifting  of  the  load  beam  on  the  dynamometer. 
The  speed  was  then  re-adjusted  and  the  spark  re-adjusted 
again  for  maximum  torque. 

The  throttle  was  run  fully  open  and  the  fuel-air 
ratio  varied  by  adjusting  the  gas  flow  rate  only. 

An  exhaust  gas  sample  was  introduced  into  the  gas 
analyzer  at  this  time,  and  the  timer  initiated  to  re¬ 
cord  the  time  and  number  of  engine  revolutions  for  the 
burning  of  0.500  cu.  ft.  of  gas. 

During  the  test  run  the  exhaust  gas  was  analyzed 
for  CO ^  f  0^  and  CO  and  the  following  data  recorded? 


Test  No. 

Oil  Pressure 
Compression  Ratio 
No.  of  Revolutions 


Cooling  Water  Temperature 
Barometric  Pressure 
Spark  Angle 
Time 


Fuel  Temperature 


Gas  Meter  Pressure 
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Air  Temperature  Pressure  of  Air  Manometer 

Cylinder  Temperature  %  C00  by  Volume 

%  by  Volume  %  CO  by  Volume 

Dynamometer  Load 

The  above  was  repeated  for  a  series  of  "runs"  at  the 
same  compression  ratio  at  different  fuel  flows ,  Immediate¬ 
ly  after  a  series  of  "runs"  the  engine  was  "motored"  to 
determine  the  friction  horsepower,,  A  typical  set  of  data 
is  shown  in  Appendix  CD 


3,2-2  Calculations 


1,  In  order  to  compensate  for  variations  in  ambient 
conditions  a  correction  was  made  to  refer  all  results 
affected  by  ambient  conditions  to  the  standard  conditions 
adopted  by  S»A,E,  i,e,  atmospheric  pressure  of  29,00  in, 

Hg.  and  85  °F,  The  correction  factor  is  calculated  by 

(3,2-1) 


the  following  formula  CF  = 


29,00 


460  +  T 


Bar,  Pres 


460  +  85 


15 


2,  Flow  of  Fuel 

To  determine  the  mass  rate  of  flow  of  propane  the 
following  equations  were  used? 


Gas  Pressure 


29,92 

ins,  of  H20  x  33o89g”x  12 


(ins,  of  Hg,) 


(3,2-2) 
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Saturation  Pressure  of  H^O: 


sat 


oq  qo 

Sato  Pres .  +  (i^s0  Hg)  (3„2-3) 

14o696 


where  Sato  Pres „  was  obtained  from  Steam  Tables 


16 


Total  Pressure  (P  )  =  Bar.  Presc  +  Gas  Pres«  (ins0  Hg) 


T 


Volume  (V)  ”  V 


st 


o  P 


st 


T 


from 


(Pf  -  P  J  (460  +  60)  , 

Vo  29.92  (460  +  Gas  Temp.)  (ft  J  <3'2~4) 


Using  the  specific  volume  of  propane  at  60  F  and  1 


atmosphere  =  8.,  503  ft  /lb 


18 


it  is  possible  to  calculate 


the  weight  of  propane  flowing  per  hour. 


V  x  60 

Therefore,  lbs.  of  CjHg  per  hour  =  gTsoTim^ 


( lbs/hr ) 


(3.2-5) 


3  o  Torque 

|  n  f ' 

The  torque  developed  by  the  engine  was  obtained  from 
the  product  of  the  load  on  the  dynamometer  scales  and  the 
moment  arm  of  the  dynamometer  (10»52  ins) 0 


10 o  52 

12 


Torque 


x  load  (ft,  lbs) 


(3  o  2  -6 
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4„  Brake  Horsepower 


BHP  = 


2ir  (Torque)  rQ  p.m. 

33,000 


(H6  P0  ) 


(3  „  2-7 ) 


5o  Friction  Horsepower 


FHP 


(Fro  Load)  r<,  Pom, 

!  33,000 


(H0  P. )  (3  o  2-8) 


6o  Indicated  Horsepower 

IHP  =  BHP  +  FHP  (HoPo) 

7o  Corrected  Indicated  Horsepower 
CIHP  -  CF  o  IHP  (Ho  Pc) 

8o  Corrected  Brake  Horsepower 

CBHP  =  CIHP  -  FHP  (HoPo) 

9o  Thermal  Efficiency 

The  thermal  efficiency  is  the  ratio  of  the  work  out 
to  the  heat  added,  expressed  as  a  percent «,  Since  this 
combustion  process  is  at  constant  volume  and  the  ^0  formed 
is  in  the  gaseous  state,  the  constant  volume  Lower  Heating 
Value  was  used  as  the  available  heato  LoHcV  -  19,929  BTU/lb. 


12 
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Indicated  Thermal  Efficiency  (tj)  =  ^ 

Brake  Thermal  Efficiency  (t^)  =  ib/hr!  Q9^  92  9)  ^ 


10c  Mechanical  Efficiency 


Mechc  eff 


CBHP 
CIHP  X 


100 


11 o  Fuel  per  Indicated  Horsepower 


lb,  fue 1/Hr/C IHP 


lb, 

CXHP-Hr , 


12c  Fuel  per  Brake  Horsepower 


lb.  fuel/Hr/CBHP 


lb. 

CBHP -Hr. 


13.  Indicated  Mean  Effective  Pressure 


IHP 


PLAN 

33,000 


(by  definition) 


Therefore 


IHP (33 , 000) 

LAN 


where 


L  =  length  of  stroke 
A  =  area  of  piston 

N  =  number  of  power  strokes  per  minute 

r .  Poim 

2 


In  order  to  compensate  for  ambient  conditions  multiply 


. 


-I* 


' 
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by  the  correction  factor  CF  or  rather  than  use  1HP  use 


CIHPo 

Therefore 


XMEP 


(CF) IHP (33,000) 

LAN 


14c  Fuel  to  Air  Ratio 


From  the  Ideal  Gas  Law,  PV  =  mR  T 


V 


mR  t 
P 


_P 

mR’T 


and  R  for  air  -  53,34  ft-lbf/lbm,  °R 


Air  density 

P  s 


Bar, 


Head  of  air  (h) 


Pres , 

14,696 

144 

(lb/ft3) 

29,92 

53,34  •  Air  Temp, 

(h) 

ins,  on 

manometer  (62,4) 

(ft,  of  air) 

Velocity  =  C  \j  2gh  =  0,062  ]j  64,4  h 


where  C  is  the  orifice  coefficient, 
Vel  (Area)  60  (Density) 


Flow 


Fuel  Air  Ratio 


144 


lb  fuel  per  hour 

lb  air  per  hour 


CF  (Ib/hr) 


A  sample  calculation  is  shown  in  Appendix  D. 

A  computer  program,  using  Fortran  IV  language,  was 
written  by  the  author  to  do  the  above  calculations.  Besides 
doing  the  above  calculations  the  program  was  designed  to 


•V 

, 

i'* 
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calculate  the  fuel-air  ratio  from  the  exhaust  gas  analysis , 
as  a  check  on  the  exhaust  gas  analysis e  This  program  also 
rearranges  all  results  in  descending  order  of  fuel-air 
ratio  before  printing  them,  and  it  also  plots  some  of  the 
results  against  fuel-air  ratio.  This  program  will  also 
calculate  the  similar  results,  if  the  fuel  is  gasoline. 

The  program  is  shown  in  Appendix  E. 


CHAPTER  IV 


RESULTS 


4.1  EXPERIMENTAL  RESULTS 

A  summary  of  the  more  pertinent  results  are  tabu¬ 
lated  in  Table  4.1.  A  sample  set  of  results,  as  calcul¬ 
ated  by  the  computer  program  is  shown  in  Appendix  F. 

The  results  tabulated  in  Table  4.1  are  shown  plotted 
in  FigSo  4.1  to  4„4  inclusive. 

Fig.  4.1  shows  the  corrected  (to  Standard  Atmospheric 
conditions)  Indicated  Mean  Effective  Pressure  plotted 
against  fuel-air  ratio.  As  one  would  expect  the  m.e.p, 
decreases  as  the  fuel-air  ratio  decreases.  Increasing 
the  compression  ratio  increases  the  m„e0pc 

Fig.  4.2  is  a  plot  of  Indicated  Thermal  Efficiencies 
at  various  compression  ratios  plotted  against  fuel-air 
ratio.  Increasing  the  compression  ratio  increases  the 
thermal  efficiency.  As  the  fuel-air  ratio  decreases  the 
thermal  efficiency  increases  to  a  maximum  value  and  then 


decreases . 


■ 
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Fig.  4.3  shows  the  effect  of  fuel-air  ratio  on  Brake 
Thermal  Efficiency.  Increasing  the  compression  ratio  has 
a  pronounced  effect  on  increasing  the  brake  thermal  effi- 

cx^iicy . 

Fig.  4.4  shows  "Comparative  Fuel  Consumption  Loops" 
showing  the  relationship  of  Indicated  Thermal  Efficiency 
and  Indicated  Mean  Effective  Pressure  and  the  effect  of 
compression  ratio. 


4.2  THEORETICAL  RESULTS 


Fig.  4.5  shows  a  comparison  of  theoretical  thermal 
efficiencies  plotted  against  percent  theoretical  air  for 

the  results  as  obtained  in  this  work  with  the  earlier 

r 


work  of  Goodenough  and  Baker 


9 


1 


,  and  results  from  Taylor 


which  he  obtained  by  using  the  charts  of  Hottel  ,  Williams , 


and  Satterfield.  The  latter  two  sets  of  results  are  for 
gasoline  as  the  fuel  and  are  not  primarily  interested  in 
very  lean  mixtures  but  range  from  lean  to  rich. 

Figs.  4.6  and  4.7  show  the  effects  of  compression 
ratio  and  fuel-air  ratio  on  the  Theoretical  Mean  Effective 
Pressure  and  the  Adiabatic  Flame  Temperature  (T^)  ,  respectively, 


'  fc  J  • 


INDICATED  MEAN  EFFECTIVE  PRESSURE  (psi) 


40 


FUEL  AIR  RATIO  (LB  «>  FUEL/LB .  AIR) 


FIGo  4.1  INDICATED  MEAN  EFFECTIVE  PRESSURE 


INDICATED  THERMAL  EFFICIENCY  (%) 
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FUEL  AIR  RATIO  (LB.  FUEL/LB.  AIR) 


FIG s  4.2  INDICATED  THERMAL  EFFICIENCIES 


BRAKE  THERMAL  EFFICIENCY  (% 
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.03  .04  .05  .06 

FUEL  AIR  RATIO  (LB.  FUEL/LB,.  AIR) 


FIG.  4.3  BRAKE  THERMAL  EFFICIENCY 


* 


INDICATED  THERMAL  EFFICIENCY  (%) 
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INDICATED  MEAN  EFFECTIVE  PRESSURE  (psi) 


FIG.  4.4  "COMPARATIVE  FUEL  CONSUMPTION  LOOPS" 


THEORETICAL  THERMAL  EFFICIENCY  {%) 
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PERCENT  AIR 


FIG .  4 „ 5  THEORETICAL  THERMAL  EFFICIENCY 
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MEAN  EFFECTIVE  PRESSURE  (psi) 
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FIG.  4.7  ADIABATIC  FLAME  TEMPERATURE 
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CHAPTER  V 


DISCUSSION  OF  RESULTS 


5.1  GENERAL 

These  experiments ,  as  previously  mentioned,  show  con¬ 
clusively  that  propane  is  an  ideal  fuel  for  use  in  high 
compression,  spark  ignition,  internal  combustion  engines. 

The  primary  object  of  the  experiments  was  to  obtain 
a  relationship  between  thermal  efficiency,  compression 
ratio  and  fuel-air  mixture  strength  for  lean  mixtures  of 
propane  gas  as  a  fuel  in  a  spark  ignition  internal  com¬ 
bustion  engine.  With  reference  to  Fig.  4.3  it  is  quite 
obvious  that  maximum  brake  thermal  efficiency  regardless 
of  compression  ratio,  occurs  at  approximately  0.04  lb.  of 
fuel  per  lb.  of  air,  i„e.  at  160%  theoretical  air,  for 
this  particular  engine  running  at  60  G;  r.p.m. 

Fig.  4.2  shows  that  increasing  the  compression  ratio 
has  a  smaller  effect  on  increasing  the  indicated  thermal 
efficiency  than  on  increasing  the  brake  thermal  efficiency. 
The  maximum  thermal  efficiency  occurs  at  a  leaner  mixture 
than  does  the  brake  thermal  efficiency.  Increasing  the 
compression  ratio  causes  the  maximum  indicated  thermal 
efficiency  to  occur  at  a  leaner  mixture  due  to  the  increase 
in  the  fuel  reaction  rate  as  a  result  of  the  increased 


pressure . 
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Fig.  404  indicates  as  one  would  expect  that  the 
maximum  thermal  efficiency  does  not  occur  at  maximum  mean 
effective  pressure  but  at  approximately  100  to  110  psi 
regardless  of  the  compression  ratio. 

With  reference  to  Fig.  4.4  it  may  be  concluded  that 
for  this  engine  running  at  600  r.p.m.  the  maximum  indicated 
mean  effective  pressure  of  the  engine  when  the  fuel  is  a 
lean  mixture  of  propane  at  a  compression  of  9  or  10  is  the 
same  as  the  m.e.p.  when  the  fuel  is  gasoline  at  a  com¬ 
pression  ratio  of  7.75  (The  highest  CR  at  which  the  engine 
will  operate  without  knock).  However ,  the  thermal  effici¬ 
ency  of  the  engine  when  the  fuel  is  propane  is  greater 
than  when  the  fuel  is  gasoline  at  the  lower  compression 
ratio.  Lowering  the  compression  ratio  results  in  a  re¬ 
duction  in  both  the  thermal  efficiency  and  the  indicated 
mean  effective  pressure. 

5  o  2  COMPARISON  -  EXPERIMENTAL  AND  THEORETICAL  RESULTS 

Although  the  air  standard  Otto  cycle  indicates  that 
the  thermal  efficiency  for  air  increases  to  100  percent 
(equation  2.3)  this  is  not  the  case  with  fuel-air  mixtures 
even  though  the  mixture  may  be  made  leaner,  i.e.  approach¬ 
ing  100%  air.  The  efficiency  reaches  a  maximum  value  and 
then  decreases ,  prior  to  the  leanest  mixture  on  which  the 
engine  will  continue  to  run  smoothly.  This  may  be 


■ 
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attributable  to  the  fact  that  the'  engine  must  rely  on 
the  fuel,  only,  for  its  heat,  which  is  not  the  case  when 
using  the  ideal  cycle  for  air,  An  optimum  fuel- air 
ratio  is  reached  beyond  which,  due  to  the  great  excess  of 
air,  the  fuel  does  not  ignite « 

Although  the  thermal  efficiencies  and  brake  mean 
effective  pressures  obtained  experimentally  are  lower 
than  the  theoretical  results,  due  to  losses  previously 
mentioned  in  Chapter  II,  when  the  curves  are  compared 
they  both  have  the  same  general  trends  and  experimental 
results  are  approximately  70%  of  the  theoretical  result s» 

5  o  3  EXPERIMENTAL  ERROR 

Although  as  many  variables  as  possible  were  con¬ 
trolled  and  a  correction  factor  applied  to  relate  all 
calculations  affected  by  ambient  conditions  to  a  standard 
base,  variations  in  ambient  conditions  did  affect  results 0 
As  an  example  ,  the  rate  of  heat  transfer  from  the  engine 
and  from  the  room  into  the  intake  manifold  vary  as  the 
room  temperature  changes  and  these  variations  ,  which 
could  not  be  determined  ,  did  have  an  effect  on  the  opera¬ 
tion  of  the  engine o 

Speed  fluctuated  during  tests  ,  even  though  the  set¬ 
tings  were  unaltered,.  These  flucuations  ,  at  very  lean 
mixtures  were  in  the  order  of  5%  and  although  they  would 
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have  a  minimal  effect  on  the  averaged  rates  of  fuel  and 
air  consumption  they  would  have  an  effect  on  the  torque . 

In  an  effort  to  compensate  for  this,  all  readings  were 
taken  in  the  same  order  for  each  testo  It  may  also  be 
due  to  these  flucuations  that  some  of  the  exhaust  gas 
analyses  were  incorrect . 

5.4  MISCELLANEOUS  OBSERVATIONS 

The  engine  operates  smoothly  on  fuel-air  ratios 
down  to  approximately  0.028  lb.  propane/ lb.  air  (230%  air) 
at  which  point  the  engine  begins  to  misfire  due  to  in¬ 
complete  combustion  of  the  fuel. 

The  mean  temperature  of  the  combustion  chamber 

O 

is  less  than  1000  F  in  the  range  tested,  ranging  from 

O 

900  at  the  higher  fuel -air  ratio  down  to  approximately 
500"  at  the  very  lean  mixtures.  The  mean  temperature  is 
lower  when  using  propane  than  when  using  gasoline  as  the 
fuel . 

The  mild  steel  exhaust  valve  functioned  well  during 
the  23  hrso  it  was  installed.  There  was  no  loss  of  com¬ 
pression  due  to  corrosion  and  observation  indicates  no 
apparent  corrosion  of  the  valve.  An  attempt  was  made 
to  measure  the  exhaust  valve  temperature  by  means  of  a 
thermocouple  mounted  flush  in  the  face  of  the  valve  and 
the  wires  extending  through  a  hole  drilled  axially  along 
the  valve.  In  order  to  avoid  modification  of  the  rocker 
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arm  mechanism  a  radial  hole  was  drilled  at  the  top  of 
the  valve  stem  just  below  the  rocker  arm  to  bring  out  the 
thermocouple  wires.  This  procedure  did  not  work  due  to 
the  fact  the  valve  rotates  and  the  thermocouple  wire  in¬ 
sulation  frayed  and  shorted  the  wires. 
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CHAPTER  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  CONCLUSIONS 

This  thesis  has  described  a  theoretical  and  experi¬ 
mental  investigation  of  the  burning  of  lean  mixtures  of 
propane  gas  in  a  spark  ignition  internal  combustion  engine 
It  has  been  shown  that  lean  mixtures  of  propane  are  satis¬ 
factory  for  the  operation  of  spark  ignition  internal  com¬ 
bustion  engines. 

Although  the  theoretical  analysis  indicates  higher 
thermal  efficiencies  and  higher  mean  effective  pressures 
than  obtained  experimentally,  both  sets  of  results,  in 
general,  have  the  same  trends.  The  reason  that  the  re¬ 
sults  obtained  from  the  theoretical  analysis  are  higher 
than  the  results  obtained  experimentally  may  be  attri¬ 
buted  to  the  fact  that  the  theoretical  analysis  is  for  an 
ideal  fuel-air  cycle  i.e.  the  processes  are  assumed  to 
be  adiabatic  with  isentropic  compression  and  expansion 
and  the  effects  of  pressure  losses  during  the  induction 
and  exhaust  processes  are  ignored.  The  mean  effective 
pressures  obtained  experimentally  are  approximately  70% 
of  the  values  obtained  by  the  theoretical  analysis.  The 
theoretical  thermal  efficiency  shows  a  continuous  increase 
reaching  a  maximum  of  48.5%  at  a  CR  =  11  and  200%  air 
(The  limit  of  this  calculation) .  The  rate  of  increase 
decreases  as  the  mixture  becomes  leaner.  The  indicated 
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efficiency,  obtained  experimentally,  initially  shows  an 
increase  as  the  air  is  increased  but  it  reaches  a  maxi¬ 
mum  value  of  40o5%  at  CR  ~  gc8  and  then  decreases  quite 
rapidly,,  The  fact  that  the  thermal  efficiency  reaches 
a  maximum  and  then  decreases  may  be  attributed  to  the 
effect  of  heat  transfer  from  the  combustion  chamber,,  As 
the  fuel  air  mixture  becomes  leaner  a  fuel-air  ratio  is 
reached  where  the  amount  of  heat  transfer  becomes  greater 
than  the  amount  of  heat  available  to  do  work  on  the  pis¬ 
ton  thus  causing  a  decrease  in  the  thermal  efficiency., 
Increasing  the  compression  ratio  causes  an  increase  in 
both  the  thermal  efficiency  and  the  mean  effective  pres¬ 
sure  in  both  the  theoretical  analysis  and  experimental 
analysis o  Increasing  the  compression  ratio  also  causes 
the  maximum  thermal  efficiency,  as  obtained  experimentally, 
to  occur  at  a  leaner  mixture,  due  to  the  combustion  pro¬ 
cess  approaching  a  constant  volume  burning,  in  a  shorter 
time  interval,  thus  reducing  the  amount  of  heat  transfer . 

Maximum  fuel  economy,  as  indicated  by  the  brake 
thermal  efficiency,  is  obtained  when  the  fuel-air  ratio 
is  at  Go  04  Ibo  of  propane  per  lb„  of  air  (160%  of  theor¬ 
etical  air)  at  an  engine  speed  of  600  r,p„m.  The  com¬ 
pression  ratio  has  a  negligible  effect  on  the  maximum 
brake  thermal  efficiency  and  maximum  fuel  economy. 

Lean  mixtures  of  propane  and  gasoline,  both  at  the 
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highest  compression  ratio  without  knock,  result  in  approxi¬ 
mately  the  same  maximum  mean  effective  pressure, 132  psi. 
However,  since  propane  is  knock  free  at  higher  CR 8 s  the 
thermal  efficiency  of  propane  at  a  higher  CR  is  greater 
than  gasoline  at  the  highest  CR  it  may  be  used  (7.75). 

Lean  mixtures  of  propane  result  in  a  lower  engine 
operating  temperature  than  when  the  engine  is  burning 
gasoline. 

There  is  no  apparent  oxidation  of  the  exhaust  valve 
when  the  engine  is  operating  on  lean  mixtures  of  propane „ 


6.2  APPLICATIONS 

When  designing  an  engine,  or  increasing  the  CR  of 
any  existing  engine,  a  considerable  saving  in  operating 
costs  may  be  realized  by  using  lean  fuel-air  mixtures 
of  propane  as  the  fuel  which  also  results  in  fewer  oil 
changes  and  less  engine  wear  than  could  be  obtained  if 
using  gasoline  as  the  fuel,  in  a  lower  CR  engine.  If  the 
engine  is  operated  at  maximum  thermal  efficiency  for  a 
fixed  speed,  less  fuel  will  be  required  thus  resulting 
in  a  further  saving.  At  times  when  greater  power  is 
required  i.e.  for  accelerating  there  will  be  a  momentary 
reduction  in  thermal  efficiency  but  the  fuel-air  ratio 
may  then  be  re -adjusted  to  again  give  maximum  thermal 


efficiency. 
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6 o  3  RECOMMENDATIONS 

Further  work  should  be  done  in  this  area  to  determine 
the  effect  of  engine  speed  on  the  maximum  thermal  efficiency . 

Comparative  tests  at  part  throttle  should  also  be 
investigated. 

Since  inlet  mixture  temperature  has  an  effect  on 
thermal  efficiency  a  wortnwhile  increase  in  thermal  effici¬ 
ency  may  be  realized  by  using  the  latent  heat  of  the 
liquified  propane  to  cool  the  inlet  fuel-air  mixture. 

Since  engines  using  lean  mixtures  are  difficult  to 
start  in  cold  weather  an  investigation  should  also  be 
conducted  to  determine  the  relationship  between  the  fuel- 
air  ratio  and  minimum  fuel-air  temperature  for  starting, 
and  for  very  cold  weather  operation. 

In  order  to  avoid  misfiring  at  very  lean  mixtures, 
and  to  see  if  the  maximum  thermal  efficiency  may  be  fur¬ 
ther  increased  at  leaner  mixtures,  the  effect  of  a  strati¬ 
fied  charge,  i.e.  rich  at  the  spark  plug  and  lean  through¬ 
out  the  combustion  chamber,  should  be  further  investigated 
using  propane  as  the  fuel. 

The  theoretical  analysis  should  be  further  improved 
to  take  into  consideration  the  effects  of  heat  transfer  be¬ 
tween  the  cylinder  wall  and  charge  during  the  entire  cycle 
rather  than  assume  an  adiabatic  process.  Consideration 
should  also  be  given  to  the  effect  of  the  pressure  losses 
involved  during  the  induction  and  exhaust  processes. 
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APPENDIX  A 

SAMPLE  THEORETICAL  CALCULATION 


Based  on  the  equations  of  Chapter  II  and  assuming 
100%  theoretical  air,  a  compression  ratio  of  9=0  and  the 
fuel  and  air  entering  at  536=6 °R0 


Note 


s  and  <t>  9  s  are  taken  from  table 


11 


if  listed 


in  the  table D 


For  the  initial  calculation  assume  the  residual 


gas  fraction 


027  and  T 


5 


2300  Ro 


from  equation  2=2-2 


C3Hg  Hr  502  +  18o8N2  +  =  027  (3CC>2  +  4H2  +  18„8N2) 


3(1=027)C02  +  4(1=027)H20  +  18= 8 (1. 027)N2 

Moles  of  charge  =  1  +  5  Hh  18=8  +  0=27(3  +  4  +  .18=8) 

=  24=8  +  =027 (25=8) 

=  24=8  +  =697  =  25=497  moles 

Moles  of  products  -  (3  +  4  +  18=8)1=027 

=  (25=8)1=027  =  26=497  moles 


from  equation  2=2-3 


U. 


.  +  H  +  HL-  +  U  -  M(1=986)T1  1  -  — 

fuel  02  N2  res „ gas  1  \  CR, 


5(83)  +  18,8(81)  +  23=8(1066)  +  180  +  1066  +  U 


res  =  gas 


. 


. 


a 


59 


-  M  (1 . 986)  T][ 

28,551  +  U  -  M (1. 986) T.  (  f  ) 
res  1  \^9  J 


Energy  of  residual  gas 


CO  2 

=  3 ( . 027 ) 

17,680 

=  1430 

H2° 

=  4 ( . 027 ) 

12,835 

«  1388 

% 

«  18. 8 ( . 027 )  9,817 

~  4980 

U 

=  7798  BTU 

res 


Try  T^  -  600 

tL  =  28,551  -f  7798  -  25.497  (1.986)  §(600) 

1  9 

=  36,349  -  45.0(600) 

=  36,349  -  27,000  =  9,349  BTU 


Considering  mixture  at  600 


U  =  o  272  (600  -  520)  +  .0160(6002  -  5202)  -  - 

C3W8  3 

(600 2  -  5202) 


272(80)  +  .0160(1120)80  - 


543 


3  x  10“ 


(216 


21.8  +  143.3  -  134.0 


543 
x  10^ 


-  141) 


1412.8  BTU 


60 


u 

°2 

-  5  (402) 

2010  =  2010 

CM 

p2 

=  18,8(1, 027 ) 395  = 

(19,3)395  =  7430 

°co2 

=  3  (.  027)  570 

(.081)570  =  46 

uh2o 

=  4(, 027)490 

— 

(.108)490  =  53 

)  U,  =  10,952  BTU 

Try  = 

580  ° 

1 

ui 

=  36,349  -  45, 

0(580)  * 

=  36,349  -  26,100  =  10,249  BTU 

Consider 

mixture  at  580° 

"3H8 

=  ,272(60)  +  , 

543 

0160(1100)60  -  — 2^4™.  (194  _  i4i) 

3  x  10" 

=  16,3  +  105,9 

-  97,7 

=  977.6 

U 

°2 

=  5(301) 

1505 

=  19,3(295)  = 

5690 

°«2 

=  ,081(424)  = 

34,4 

V 

=  ,108(490)  = 

53 

U  =  8260,0  BTU 


o  T^  is  between  580  and  600  °R, 
tion  to  get  T^. 


Do  a  linear  interpola- 


•  v'*' 


61 


DT 


20  -  DT 


10,249  -  8260  10,952  -  9349 


DT 


20(1989) 

1989  +  1603 


20(1989) 

3592 


llo  1 


T  =  580  +  11  =  591  R 


Adiabatic  compression  from  T^  to  T2 


$ 


C3H8 


rqi 

272  In  +  o  0320  (591  -  520! 

o  z  u 


543  (5912  -  5202) 

.  o 


2  x  10 


272  In  1,138  +  ,0320(71) 


543 


2  x  10 


8 


(1111) 71 


0351  +  2,27  -  ,214 


2,091 


O 


2 


5(  ,56  +  (.75  -  56; 


5 (. 56  +  ,10)  =  5  (.66)  =  3,30 


<P 


N, 


19, 3 ( o  54  +  —-(.72 


54)  =  19,3 (, 54  +  ,10) 


19 , 3  ( „ 64)  =  12,33 


<t> 


CO 


2 


08l(  .77  +  |~(1,Q2 


77)  )  =  . 081 ( , 77  +  .138) 


081 (.908) 


07 


4> 


H2° 


108^.71  +  “(.94  -  »71) 


108  (.71  +  .126) 


108 ( . 836)  =  .09 


<P 


591 


17 . 88 


62 


Try  =  1200 


<D 


CnH 


8 


.  272  In  ^  +  „  0320  (1200  -  520) 


543 


3  x  10 


8 


(12002  -  52Q2 


272  In  2-31  +  ,0320(680)  - 


543 


3  x  10 


8 


(1720)680 


,228  +  21,8  -  3,18 


*=• 

18.85 

■e- 

O 

ii 

5(4,50) 

22.50 

<t> 

N2 

19.3(4.25)  = 

0 

0 

0 

ro 

00 

C02 

.081(6.91) 

.56 

H«0 

.108(5.31) 

.57 

z 

y  <i>i2oo = 

125.48 

As  -  1 1?00  -  Y  *591  -  1.986(H)  in  CR 

=  125,48  -  17,88  -  1,986(25,497)  In  9 

=  125,48  -  17.88  -  111,2 

=  125.48  -  129.08  =  -3,60  BTU 

. assumed  too  low 


Try  T2  ™  1300 


H 

3  8 


272  In  +  .0320  (1300  -  520) 


543  2  2 

(1300  -  520  ) 


2  x  10 


8 


520 


63 


=  .272  In  2 „ 5  +  .0320(780) 


543 

2  x  108 


. 249  +  24.9  -  3.86  =  21.29 


4)  =  5(4.98) 

U2 

<t>N  =  19.3(4.68) 

4)co  -  .081(7.69) 

(DH  =  .108(5.85) 


24.9 

90.2 

.62 

.63 


1300 


137.64 


As  *  137.64  -  17.88  -  111.2 

=  137.64  -  129.08  =  8.56  BTU 


Linear  interpolation  to  get  TL 


DT 


100  (3.60) 

3.60  +  8.56 


360 

12.16 


29.6 


=  1200  +  30  *  1230 °R 

Determine  the  energy  of  mixture  at  1230° 

2  543 

U  =  .2721'  +  .0161  -  — 1 

38  3  x  10 

=  .272(1230  -  520)  +  ,016(12302  -  520 


(12303  -  5203) 

.272(710)  +  .016(1750)710 


543 

3  x  102 


(1820)780 


543 

3  x  108 


(1860  -  163) 


\  • 


64 


U 


O . 


U. 


N, 


U 


CO, 


U. 


H2° 


CE 


C  H 
3  8 


.  272  (710)  +  .016 (1750)710  - 
193  +  19,870  -  3077  = 


543 


300 
16,986  BTI.J 


(1697) 


5^3675  +  . 3 (4263  -  3675) 


5 (3851) 


19,240 


V  ufc  I  f  W  > 

.  108(^4339  + 

878,822  B/mole 


3 (3996 

-  3461)  J 

=  19. 3 (3622) 

=  69,900 

3  (6721 

-  57  36) 

=  .081(6032) 

=  488 

3  (5030 

-  4339)  ) 

=  .108(4546) 

=  491 

U2  =  107 , 105BTU 


Total  energy  =  878,822  +  107,105 


985,927  BTC 


Substitute  into  equation  2.2-11  and  add  like  terms 


Assume  T^  ~  5200  R 


Product 
3. 081 (x) CO 
3.081 (1-x) CO 2 
4. 108 (y) H2 
4. 108(l-y)H2O 


Numerical 

Terms 


166,000 


172,500 


x  terms  y  terms 
90,200 
-166,000 

111,200 

-172,500 


|-(x)  +  |-(y) 


(1.027)02 


47,900  63,900 


19. 3N2 


559,000 

897,500  -27 , 900 (x)  2,6QQ(y) 


o  a,  ;rj 

I ! 
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CEco  = 

-  3 (1.027) (x) 

121,181 

=  373,000 

(x) 

CV  ■ 

=  4(1.027) (y) 

103 ,486 

=  424,500 

(y) 

Energy  balance  gives?  - 


985,927 
.  88,427 


897,500  -  27, 900 (x)  +  2,600(y)  +  373,000 (x)  +  424,500 (y) 
345 , 100 (x)  +  427  100 (y) 

88,427  -  345 , 100 (x) 

427,100 


0 o  2065  -  Q.808x 


From  the  tables 


12 

_ 


lo9  5*3  =  °'829 


=  6.75 

Substitute  into  equation  2.2-13 


y  * 


_ _ x _ _ 

6.75  -  5.75x 


Equating  y's 

y  -  0.2065  -  0.808x 


__ _ x _ _ 

6.75  -  5 . 75x 


©  e 


(0.2065  -  0 . 8Q8x)  (6.75  -  5.75x)  =  x 

1.393  -  5.45x  -  1.189x  +  4.65x2  =  x 


4 . 65x  -  7.639X  +  1.393  =  0 


x 


7.639  -  V  (7.639)^  -  4(4.65)1.393 

2 (4.65) 


*:  : 


7.639 


66 


-  \l 58.3  -  25 „  9 

2(4.65) 


7.639  -  5.69 

2  (4.65) 


1.949 

2(4.65) 


0.2097 


y  =  0.2065 

=  0.2065 

Obtain  P0  from  P0 


0.808(0.2097) 
0.1692  =  0.0373 

cr(t2) (p1) 


9(1230) (1) 
2  591 


18.73  atm. 


Obtain  KCQ  ,  using  equation  2.2-8 
2 


Kc°2 


|(l+f)x  +  |(l+f)y 


/  0.2097  V 

\1  -  0.2097 ) 


1.0  81 
2 


(.2097) 


+ 


4.108  t 
2  { 


18.73 (5200) 

25.497 (1230) 


0.2097V 
0.7903 ) 


0.322  +  0.0765 


18.73 (5200) 

25.497  (1230) 


=  (.265) 2 


(.3958) 


18.73 (5200) 

25.497 (1230) 


.0861 


From  the  Tables 


log  K 


=  -1.32 


0373) 


.  ■  ..A1 
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KCo2  =  _0^_048 


assumed  too  low 


Try  =  5400 


Product 


Numerical 

Terms 


3(1.027) (x) CO 
3 (1.027) (l-x)CO, 


174 , 000 


4(1.027)  (y) H 


2 


4(1.027)  (l-y)H20 


182,000 


i<*>  +  f 


(1.027)  O, 


19.3  N 


2 


585,000 


Energy  balance  gives: 


x  terms  y  terms 


94,300 
•174, OOP 


117,000 

•182,000 


50,200 


66,900 


941,000  -29 , 500 (x)  l,900(y) 


985,927  =  941,000  -  29,500(x)  +  l,900(y)  +  373,0Q0(x)  +  424,500(y) 

.  44,927  =  343 , 500 (x)  +  426,400(y) 


Y 


44,927  -  343, 500 (x) 

426,400 


=  0.1052  -  0 . 805 (x) 


From  Tables 


12 


log  =  0.840 


XWG  “  7*06 


68 


.  .  y 


_ x _ 

7.06  -  6  o  06 (x) 


Equating  y's 


y 


0  . 1052 


0 . 805x  = 


_ x _ 

7 . 06  -  6 . 06x 


(0.1052  -  0 o 805x)  (7.06  -  6.06x)  =  x 

0.743  -  5.68x  -  0.638x  +  4.88x2  =  x 
4 o 88x2  -  7.318x  +  0.743  =  0 

7.318  -  1/(7.318)2  -  4(4.88)0.743 

x  ~  . 

2  (4.88) 

=  7.318  -  1/53.5  -  14.5 

2(4.88) 


7.318  -  6.24 

2  (4.88) 


1.078 

2 (4.88) 


0. 1102 


y  =  0.1052  -  0.805(0.1102) 


0.1052  -  0.0887  =  0.0165 


|81-  (.1102)  +  (.0165) 


18.73 (5400) 

25.497 (1230) 


K 


0.1102 


CO. 


0.1102 


/0.1102\2 

\^0 . 8898 J 


0.1699  +  0.0338 


18.73 (5400) 

25.497 (1230) 


=  (0 . 125) 2 


0.2037 


18.73 (5400) 

25.497 (1230) 


=  0.01027 


69 


From  Tables 


12 


lo9  KC0, 


■0. 94 


K  =  0=115 

u2 


T. 


5400  R  is  too  high 


Do 


a  linear  interpolation  to  obtain  and 


DT 


200  -  DT 


0=0861  -  0=048 


0=115  -  0=01027 


DT 


0=0381 (200) 

0=0381  +  0=10473 


0=0381 (200) 

0=14283 


=  53  =  4 


T. 


5200  +  53=4 


5253=4  R 


x. 


0=2097  - 
0=1831 
0=0373  - 


53  =  4 
200 


53  =  4 
200 


(0=2097  -  0.1102) 


2097  -  =0266 


=0373  -  0=0165) 


0373  -  =0055 


0=0318 


Entropy  at  T^ 


Using  equation  2=2-14 


4) 


M  CE 

nr  ,  ,  MCO  CEC0  ,  H2  H2 

MR  In  pp  +  - - -  + - - - 


0 


CO, 


3=081(1  -  =1831) 


2^3 

27=98  +  ~r  (29=  51  -  27.98) 
500 


3=081 (0=8169) 


27.98  +  0=774 


. 
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3.081 (0=8169)  (28.754)  =  72.2 


0 


CO 


3.081 (0.1831) 


3.081 (0.1831) 


2  51 

17 . 72  +  |~~(18.57  -  17.72) 


17.72  +  0.43 


H2° 


<t> 


H, 


=  3.081(0.1831) (18.15)  =  10.22 


4.108(1  -  .0318) 


?  S  ? 

22.42  +  ~~  (23.65  -  22.42 
500 


4.108(1  -  .0318)  (22.42  +  0.622) 

4.108(0.9682)  (23.042)  =  91.6 


4.108(0.0318) 


2  S  ? 

16.81  +  (17.65  -  16.81) 


4.108(0.0318)  (16.81  +  0.424) 

4. 108(0.0318) (17. 234)  =  2.25 


0 


O. 


(0.1831)  +  ---^-(0.0318) 


2  R 1 

18.49  +  -—-(19. 38  -  18.49) 


(0.282  +  0.0645) (18.49  +  0.45) 
0.3465(18.94)  =  6.57 


0 


N, 


=  19. 3 


2  6? 

17.54  +  (18.38  -  17.54) 

500 


19.3  (17.54  -  0.424)  =  19.3(17.964)  -  346.8 


0  =  529.64  BTU 


M. 


T. 


M 


2 


T. 


I) 801 > 


'  M.t  '  '5  ' 


71 


M3  =  3  o 081  +  4  o 108  +  19.3  +  0.3465  =  26.8355  moles 


=  26.8355 

3  25.497 


5253 

1230 


18.73  =  84.3  atm. 


Constituent 

££ 

In  pp 

1.986(M)ln  pp 

3.081 (0. 8169)C02 

7.90 

2.065 

10.31 

3.081 (0  = 1831) CO 

1.77 

0.572 

0.64 

4.108(0. 96 82) H20 

12.46 

2.52 

19.83 

4 . 108 (0 . 0318) H2 

0.41 

-0.891 

-0.231 

0, 346502 

1.09 

0.0862 

0.059 

19. 30N2 

60.60 

4.11 

157.4 

I 

188.01 

s  3  “  529.64 

-  188.01 

+ 

3. 081 (. 

1831)121,181 

5253 

4 . 108  ( . 0318) 103 

,486 

r 

5253 

=  529.64 

-  188.01 

+ 

13.00  + 

2.57 

=  357.20 

BTU 

To  find 

Assume  =  3300  °R 

from  Tables , 

kC02  = 

12  ,  log 

8.32  x  10 

KCO  =  _7-08  and  109 

'8  =  3*92 

.593 


72 


Try  x  =  0 . 006 


Substituting  in  equation  2 . 2-15 


K 


CO. 


—  /  Q  ° ooo 

0.994 


2  r 


3 „ 081 ( . 006) 


4. 108  (.006) 

2  (3 . 92  -  2.92  (.006) ) 


3300 


25 . 497 (591) 


=  4.  01  x  10 


-5 


0.00924  + 


4.108(0.003) 

3.92  -  0.0175 


3300 


25.497 (591) 


4.01  x  10 


-5 


0.00924  +  0.00315 


0.219 


4.01  x  10' 


0.01239 


0.219 


1.09  x  10 


Try  x  -  0.0057 


Kco2 


QaQ057\2 

0.9943/ 


3.081  (.0057)  4  0 108  ( . 0057 ) 

2  2(3.92  -  2 . 92 ( . 0057 ) 


0.219 


=  3.28 

X 

10' 

-5 

VO 

r- 

00 

o 

o 

0 

o 

L _ 1 

.  4.108(. 

0057) 

0.219 

2(3.92  - 

.0166) 

=  3.28 

X 

10' 

-5 

0.00876 

+  0.00299 

0.219 

00 

CM 

O 

ro 

1! 

X 

10' 

-5 

0.01175 

0.219 

=  8.43 

X 

10 

-8 

Try  x  =  0.00568 


0  o  00568\2 
0.99422  J 


3.081 

2 


(.00568) 


4  o 108 (. 00568) 

2 (3.903) 


0.219 


3.26  x  10 


5 


0.00874  +  0.00298 


0.219 
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3  .  26  x  10  5  (0,01172)  0,219  =  8,35  x  10  8 


so  x  =  0,00568 


0,00568 


0,00568 


3,92  -  2 , 92 ( . 00568)  3,92  -  ,0166 


0,00568 

3,9034 


0.00145 


0 


CO, 


3,081(0,99432) 


43.76 


67,1 


$ 


CO 


3,081 (0. 00568) 


28.15 


0,246 


<D 


H2° 


=  4.108(0,99854) 


34,49 


70.6 


4>  =  4,108(0.00145)  13.32 

ii^ 


0.079 


4>  =  0.01172  (14.69) 


O 


2 


0.172 


<j)  —  19,30 


27.87 


268.  8 


407.00  BTU 


M4  =  3.081  +  4.108  +  0.01172  +  19.3  =  26,501  moles 


26.501  3300(1) 

25.497  ’  591 


=  5.82  atm< 


74 


Constituent 

In  pp 

1.986(M)ln  pp 

3.081(0.99432) C02 

0.672 

-0.397 

-2.415 

3.081(0.00568) CO 

0.00384 

-5.56 

-0.193 

4.108(0. 99854)H20 

0.900 

-0.105 

-0.855 

4.108(0. 00145) H2 

0.00131 

-6.63 

-0,078 

0.01172  02 

0.00258 

-5.96 

-0.139 

19.30  N2 

4.24 

1.445 

55.4 

\l.986(M)ln  pp  = 

51.72  BTU 

407.00  -  51.72  +  ^  12 1-^°- 


4, 108(, 00145)103, 500 

3300 


=  407.00  -  51.72  4-  0.64  +  0.19  =  356.11  BTU 


assumed  too  low 


Assume  =  3400  R 


From  Tables 


12 


log  Kco  =  -6.60  and  log  =  0.616 


K 


CO. 


2.52  x  10  and 


4.13 


Try  x  =  0.008 


/  Q.0082\2 

3, 081 (.0082) 

4. 108 (.008) 

y0 . 9918 / 

2 

2(4.13  -  3.13 (.0082) 

25 

6.5  x  10 


0.0123  + 


0.0164 


4.13  -  0.025 


0.226 


3400  _ 

497 (591) 


75 


=  605  x  10 


-5 


0  o  0123  +  0 , 0040 


0,226 


K 


CO. 


=  6,5  x  10  5 

(0,0163)  0,226 

=  2,40  x  10“7 

0,0082 

/0,0082\2 

3 , 081 ( , 0082 )  . 

4, 108 ( , 0082 ) 

^0,9918/ 

2  + 

2(4,13  -  3 , 13 ( , 0082 ) 

0 


6,84  x  10 


-5 


0  o 0126  +  0,0041 


0,226 


6,84  x  10  5  (0,0167)  0,226  =  2,58  x  10  7 


Try  x  =  0,00814 


0,00814\2 
0,99186  ) 


3,081 (,00814)  4 , 108 ( , 00814) _ 

2  2  (4,13  +  3.13  (,00814) 


6,74  x  10 


0,0125  +  0,0041 


0,226 


=  6,74  x  10  5  (0,0166)  0,226 


2,53  x  10 


So  x  =  0,00814 


_ 0,00814 _  _  0,00814 

4,13  -  3 . 13 ( . 00814)  4,13  -  0,0255 


0 

0 

0 


CO 

H2° 


0,00814 

4,105 


—  0,00198 


3, 081 (. 99186)22, 32  = 


3, 081 (. 00814)14, 34 


4, 108 (,99802) 17,61 


4 , 108 ( „ 00198) 13,58 


68,2 
0,36 
72 , 2 

0,11 


226 


0,226 


I  ■: 
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4> 


0 


0 


2 


N, 


0 . 0166 (14.95) 
19.30(14.20) 

4> 


0  =  25 
2740  0 

415 . 12  BTU 


—  3.081  + 

4.108  +  19.3 

+  0.0166  - 

26.506  moles 

_  26.506 

P4  25.497 

3400  «i)  = 

591  '  ' 

5.97  atm. 

Constituent 

EE 

In  pp 

1.986(M)ln  pp 

3,081(0. 991 85) CO 2 

0.689 

-0.372 

-2.26 

3 . 081 (0  o  00814) CO 

0.00565 

-5.17 

-0.257 

4.108(0.99802)H20 

0.923 

-0.080 

-0.650 

4. 108(0. 00 198) H2 

0.00183 

i 

0 

U) 

o 

-0.102 

0 . 0166  02 

0.00374 

-5.59 

~0c 184 

19.30 

4.35 

1,47 

56,2 

52.75  BTU 


„,r  .  3.081  (.00814)  121,200 

415.12  -  52o75  +  - 1 - 2400  - - - “ 

.  4.108(.Q0198) 103,500 

+  _____  __ 


=  415.12  -  52,75  +  0,89  +  0.32 

assumed  too  high. 


363.58  BTU 


By  linear  interpolation  obtain  ,  x^  and  between  3300 
and  3400 °R. 


■ 
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DT 


357.20  -  356.11 


100  -  DT 


363.58  -  357.20 


DT 


1.09(100) 

1.09  +  6.38 

3300  +  14.6 


1.09(100' 

7.47 


3315  R 


14.6 


x. 


14  fi 

0.00568  +  00814 


00568) 


00568  + 


0.00604 


14  6 

0.00145  +  i~'^00198  "  °  00145) 


00145  + 


0.00153 


"5  - 


MR  In  pp 


No.  of  moles 

-  3.081  + 

4.108  +  19. 

30  =  26.489 

Try  Tj.  - 

2300 

0R 

0  ss 

C02 

3. 

081(16.81) 

=  51. 9 

(t)  — 

h2o 

4. 

108(13.22) 

=  54. 3 

d>  - 

N2 

19 

o 30 (10.92) 

=  210.7 

r — j 

L* 

=  316.9 

BTU 

Constituent 

RP_ 

In  pp 

1 . 986 (M) In  pp 

3.081C02 

0.1162 

-2.15 

-13.16 

4. 108H20 

0.1548 

-1.863 

-15.19 

19.3N2 

0.729 

-0.316 

-12.10 

V- 

*  -40.45  BTU 

00036 


00008 


-  (■ 


78 


.  °.  sc  =  316 o  9  +  40 o 45 

5 

.  Take  Tc  =  2300  °R 

V9  M  (T  ) 

f  =  v~  m^TcrTt5 

25  . 497  (591) 

26.489(9)  2300* 

Tj.  is  found  to  be  as  assumed  but  f  is  slightly  high  0.0275 

rather  than  0.027  as  assumed  so  the  calculations  should  be 

re-done  using  f  -  0.0275,  for  greater  accuracy. 

To  determine  thermal  efficiency 

Substitute  into  equation  2.2-17 

9_  ,  __  v 

U1  amix  600  20  Umix  600  Umix  580^ 

=  10,249  -  (10249  -  9349) 

=  10,249  -  (900)  =  10,249  -  405  -  9 , 844  BTU 

u2  =  107,105  BTU 


357.35  BTU 


0.0275 


To  find  u 


3 


u 


CO 


3.081(1  -  .1831) 


2 


53,963  +  (56,569  -  53,963) 


3.081 (. 8169) 


53,963  +  (2606) 


3.081(.8169) (53,963  +  695)  =  3 . 081 (. 8169) 54 , 658 


J  I 


79 


u 


CO 


u 


H2° 


u 


H, 


u 


O, 


u 


N, 


137/500 
3„ 081 („ 1831 ) 


29/287  +  (30,665  -  29,287) 


200 


3 . 081 ( . 1831)  (29,287  +  370) 

16,780 


3.081 (. 1831) 29,657 


4.108(1  -  .0318 


£  “5  A 

42  ,087  +  (44,307  -  42,087 


4.108(0.9682) (42,087  +  592)  =  4 . 108 ( . 9682 ) 42 , 67 9 

169,400 


4, 108 (0 o  0318) 


R  8  4 

27,156  +  ~~  (28,508  -  27,156) 


200 


=  4  . 108 ( o  0318)  (27,156  +  361) 


4 . 108 (.0318) 27,517 


3,580 


31,115  +  (32,622  -  31,115) 


0.3465 

0 o  3465  (31,115  +  401) 


0 . 3465 (31,516) 


10,900 


19.  30 


58  d 

28,961  +  (30,337  -  28,961) 


19.30  (28,961  +  367)  =  19.3(29,328) 

566,000 


u. 


904,160  BTU 


To  find  u 


u 


CO. 


29,750  + 


15 


(30,991  -  29,750) 


3.081(1  -  .00604) 


100 


80 


=  3  o 081 (0»  99395)  (29,750  +  186)  =  3 0 081 ( . 99395 ) 2 9 , 936 

=  92,900 

uco  =  3 o 081  ( o 00604)  (16,414  +  100)  =  3 „ 081 ( . 00604) 16 , 514 

=  308 


u 


H2° 


u 


H, 


u 


O. 


=  4.108(. 99847)  (22,066  +  149)  =  4. 108 ( . 99847 ) 22 , 215 

=  91,300 

=  4 . 108 ( . 00153 ) (14,921  +  91)  =  4 0 108 ( 0 00153 ) 15 , 012 

=  94_ 

( .  00604)  +  ( .  00153)  (17,386  +  104) 


u 


N 


2 


(0*00933  +  0.00314) (17,490)  =  0.01247(17,490) 

218 

19.30  (16,199  +  98)  =  19*30(16,297)  =  316,300 


u4  =  501,120  BTU 


(904,160  -  501,120)  -  (107,105  -  9,844)  .  -jrn 
44.094(19,929) 


(403,040  -  97,261) 

44.094(19,929) 


100 


(305,779) 100 
44.094(19,929) 


34.7% 


To  find  m»  e „ p. 

Substitute  into  equation  2.2-16 

mRTl  25.497 (1545)591 

V1  Pn  14.7(144) 


110,000  ft* 


81 


VJL  =  110,000 

CR  9 


12,200  ft 


3 


_  305,7 7 0  (778 ) _  =  305,77 9(77 8 ) 

(110,000  -  12,200)  (144)'  97,800(144)" 


169  psi . 


■ 
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APPENDIX  B 

PROGRAM  FOR  THEORETICAL  ANALYSIS 


$ JOB  525015  G.R.  POND 

STIME  30*3000 

$  I B JOB  POND 
SIBFTC  POND  NODECK 

c  **  otto  cycle  calculations  ** 

C 

DATA  FU1 , FU2/5HC8H18 *4HC3H8/ 

REAL  N2  P  *LS*MC0fMC02 , Mus * MW > MO 2 

K  I  =  1 

RE AD ( 5  »  51 )  XC 
I F ( XC.GT.4. )  GO  TO  2 
FU  =  FU2 
XC  =  3  • 

XH  =  8  • 

CR  =  8  •  0 
CEF  =  878822  • 

HEAT=44. 098*19929. 

UCG=180.0 
GO  TO  4 

2  KI  =2 

FU=  FU1 

XC=  8  • 

XH= 1 8 • 

CR  =  7  •  0 

CEF=2201618. 

HE  AT  =  236I 226. 

UCG=674.0 
4  XH2  =  XH/2 • 0 

X02= (XH2+2.0*XC) /2. 

XN2  =  3.76*'-X02 
WR I TE ( 6  »  5  2 )  FU 
6  PC=0 • 80 
TT 1=920. 

TT2=1350. 

F  =  0 . 03 
CZ  =  XH2 

WRITE(6»54)  CR 
M=  1 
N=  1 5 

DO  11  J=M,N 
L  I  =  0 

PE=100.*PC 
WR  I  T E ( 6  » I  3 )  PE 
M0S=PC*XC2 

IF ( PC. LT. 0.999)  GO  TO  8 
C 

C  IF  EXCESS  AIR  OR  100  PERCENT  AIR 
M02  =  XC 
MCO=0. 000 
MCO  2  =  XC 

IFIPC.LT. 1.001 )  GO  TO  12 


i  .  f  Jr.  .  J,  .J  i  :  ■  <>.'  .if  •-  l  .  A 


ukus  •  >l.j  at  it 

0  U  0  £  .  0  e 
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>2  3 1  Ott 


ano^i  aoiHi 
anon  oTig : f 

OTTO  ** 


*■  ?^iu1TA-j 3J3Y?  OTTu 
1 1 A  t  .  IT.  \  J  u  1 »  0  3  /  i 

*  »  r  .i  ■  •  ■  J ‘  I  »  j*  -J  ■  .  j  •  1 

I  =  IX 

0X  (  1  2  »  5 ) GAB • 
S  uT  ju  (  .♦=  .TO.0X  )  31 


sua=ui 

.eox 

•  8  =  H  X 

. 


.  3388^°=  -133 

. 

.  -  8  I  =  0DU 
A  OT  00 
S=  IX 


.BI3IoSS=  u3 
.  as  ;d£S  =  TAjH 
v.AY3=0DU 
.  ;\HX  =  SHX 
. s v ( DX*G«$+SHX ) =  SOX 
SOX >dT»€*SMX 
U3  (S5«d)3Tlfl'w 
8  .COM 
. j  S9  =  I T  T 
•J5£I*STT 


£  C  .  =  - 
S  H  X  =  S  - 


F  (  Ac»d)  j  T  Ifl>. 


i  =;•. 
<?  i  =n 


n < ri  oa 
c  =  i  j 
os*.  r=3H 
3H  (£  !*d) BTlflW 

S JX*DS*aOM 

8  u  Oj  <'#£?•  G  .  J.0H  Ml 

/■i  i  m  UojFjh  /tu  xlA  -3j3X_, 

jX*$OM 
oC  *o=-Oi*' 

; <*sooM 
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EXC= ( 2.*X02*PC-2.*M02~XH2  )  /  2. 

GO  TO  14 

IF  AIR  LES'3  THAN  100  PERCENT  C02  DISSOC  TO  CO 
8  M0S=PCfX02 
M02=M0S~XH/4. 

MC0  =  2 • 0* ( XC-M02 ) 

IF ( MCO.GE. 0.0000 )  GO  TO  10 
MCO=  0 . 000  )0 
10  MC02=M02“MC0/2. 

12  EXC=0. 000000000 

14  CONTINUE 
P02=PC*X02 
PN2=PC*XN? 

X  LS  =  1  •  J0+PO2+PN2 
MW  =  XC*12.0  +  XH**1.008 

INITIAL  CUND@S  ^0  FIND  T1  “  FUEL  MIXTURE  PLUS  CLEARANCE  GAS 

ASSUME  FUEL  TEMP  =  537. 

**  CONSTANT  PRESSURE  ADIABATIC  PROCESS  ** 

U1  =  FU5+ HCG+HGF-PV/U+G  ASSUME  0=0.0 
F  A  =  M  W / ( 32. 0*PO2  +  2  8. G*PN2 ) 

HCG  =  PO  1*83. 0  +  PN2*8 1. 0+( XLS-1 .0)*1 066.0 
HG  =  UCG+ 10  6  6.0 

CXHY+PC*( X02+XN2 ) =MC02*C02+XH2*H20+MC0*C0+XN2*PC*N2+EXC*02 

FIRST  TRIAL  ASSUMING  NO  DISSOCIATION  AND  FuLLuwING 
Tl=TTl+80. 

T2=  T  T2  +  18  0 • 

T  4  =  3600. 

T5=2500. 

X=  0  •  0 
Y=  0  •  0 
L  F  =  0 
LK  =  0 
AZ=MC02 
BZ=MCC 
DZ=EXC 
EZ  =  PN2 

16  MC02  =  AZ 
MCO  =  BZ 
XH2=CZ 
EXC  =  DZ 
PN  2  =  EZ 

SUM  OF  MOLE  3  ON  RIGHT  SIDE  OF  EQUATION  (PRODUCTS) 

Z= 1 ,+F 
Z A=MC02*Z 
Z  B  =  MCO*Z 
ZC=XH2*Z 
ZD=EXC*Z 
ZE=  PN2*Z 
MC02  =  Z  \ 

MCO=ZB 


£ 


.  .  .  , 

11  3 


[ 

SI 

3  L  T  M  0  3  AT 
SOX  *0^  =  509 
c  /lx  *03  =  S  kl3 
s/i«  +  +  .  i  ^ejx 

1 0  G  •  ' •*’*H .  +•;  .  .  I  *  3X*rtM 

jJ*  A/i  UJJ  .  U-i'i  _i  rt ,  ;  A  i  *•  Jj  J1  -  :  1  (JM3  c-UjhO  JAITItfl 

.V£a  =  '1M3T  J3u3  BMUaEA 
'*  vicJ..  ottS  3  1  T  A  d  A  I  <  A  j>Uio^33  TtfATSMuO  ** 

. 

(  3,  3*l  .RS  +  S  .££  )  \W^=A3 
L  .ddd*  (0.  I-cJX  )  4  .I8*3A3-k  .",>*£  0M«03H 

•  +;  3U=OH 

S^*  X.i  +  Sr  *33*v>kiX  +  u0*  OK.+uSH*S1  X  +  SO*S03M=  (  £1  X  +  SCX  )*33  +  YHX0 

o* .  A  t'  *  >.  i  T  A  i  3  d  £  c  I  u  ell  ut'  1 1 1  il)  -  3  A  j  A  I  X  T  T  3. 5-1 1  3 

.08^  I T  T  = I T 
•  ^  8 1 +£  T  T  =  ST 
.wi  d £  =  A T 
•  3  L  a  S  «  €  T 
•  —  X 

0  .o  =Y 
0  =  3  J 
C  =  XJ 
SOOM=SA 
oom=S8 
0X3  =  S<J 
£H3=S3 
SA  =  SOOM  di 
SO  =  u3M 
S3=SHX 
SO  =  0X3 
S3=StfH 

(  .  .  .•  .ai  v/.s )  HU1TAUU3  3^  Kilt  i  ;  ■. .  •  IJOH  1  Q  M  U  2 

3+ • X  =  S 
/  S*  S00M  =  AS 

S*  COM*  81 
S*SHX=0S 
S* 3X3=CJS 
S*£/:M  =  JS 
/  S  =  SO> 

6  S  =  C  OM 


I  T  Kl  .  '  ,  r.  .  j  / 

.  X  *  '1  = 

.  a\.  -  .vjf-'.ssc; 

(  SUM-OX  .  A  v  .  3  =  o0f' 

a.  uT  O  (  0  V.  •  w'  •  *  0  f /  ■  )  1  1 

C(  oc  #o*oom 

.  GW  3-- 

r  0  C  .  =  0  X  . 
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XH2  =  ZC 
EXC=ZD 
PN2  =  ZE 

13  XRS=XH2+MC0+EXC+PN2+MC02 

NO.  OF  MOLES  IN  IHARGE 
15  XCH=XRS-1, 

WRITE!  3*60)  MC02*MC0*XH2»EXC*PN2*XRS»XCH*ZA»AZ 
KO  =  Q 

17  CALL  ENERGY  ( T 5  * XC * UC02 * UCO , U02 * UN2 ♦ UH20 , UH2 , UFU ) 
UA=MCC2*( X*UCO+( 1.00000-X) *UC02 ) 

U  B  =  X  H  2  *  (  Y  *  U  H  2  *’*  (  1 . 0000 0“  Y  )  *UH20  ) 

UC= ( X*MC02/2.0+Y*XH2/2.+EXC ) *U02 
UD=PN2  *UN2+MCO*UCO 

P ( V1-V0  )  = U ( V1-V0 )  =RT/V1 < V1-V2 )  =RT(1“1/CR) 

FUH=F*  (  UA  +  UB  +  i  IC  +  UD  )  +HCG+HG 
R=  1  • 98  39*XCH* ( 1.0-1.0/CR) 

18  U1=FUH~R*T1 

CHECK  U1  BY  CONSIDERING  GAS  AT  T1 
USENS=F*JCL+( UAIR+UFU) 

CALL  ENERGY  ( T 1 , XC * UCQ2 * UCO , U02 , UN2 , UH2G , UH2 * UF U ) 
UA=MC02*( X*UCO+( 1.00000-X) *UC02 ) 

UB=XH2* ( Y*UH2+ ( 1.00000-Y)*UH20) 

UC= ( X*MC02/2 .0+Y*XH2/2 .+EXC ) *U02 
UD=PN2*UN2+MCO*UCO 
UE=P02*U02+PN2*UN2+UFU 
USENS=UE+ ( UA+UB+UC+UD ) *F 
PT  1  =  T1 

CALL  TEMP ( T1 *USENS*U1 , DIF ,KO) 

T  T  =  ABS  (  (PT1-TD/T1  ) 

I  F ( TT.GT.0.Q1 )  GO  TO  17 

IF ( K0.LT.2.0R.DIF.GT .0.0 )  GO  TO  18 

V2=V1/CR 

V1  =  1545.*T1*XCH/ (144. *14. 7 ) 

K0  =  0 


ADIABATIC  COMPRESSION  TO  T2  -  AT  CONSTANT  ENTROPY 

ENTROPY  OF  FUEL 
20  I F ( XC.LT.4. )  GO  TO  22 

SFU=5.93*A  LOG ( T  2 / T 1  )+0.061*(T2~Tl) 

GO  TO  24 

22  SFU=0.272*ALQa( T2/T1 )+.0320* (T2”Tl ) - ( 543.* ( T2**2 . 0~T 1 **2 . 0 ) / 2 00000 

1000.  ) 

24  A=  1  •  0/ T2 

b=i  .0/ n 

SC02  =  F*MC02*( l.-X ) *( 14. 3*A  LUG (T2/Tl)+6530.*(A-B)~705000.*( A**2.0~B 

1**2. 0)  ) 

SC  0=  F* ( MCO+MCU2*X ) *( 7.47*ALuG (T2/Tl)+3290.*( A“B ) “ 5 3 5 0 00 . * ( A**2 . 0“B 

1**2. 0) ) 


3S=SHX 
CIS  OX  J 
3S=StfS 

S  G  j  1 4- s  h  ,  h  +  ) X  J  +0  3K  +  GH X  =  G  £ X 

>!  A  il  i  1  JGI  i  •  . 

,  I-  '  =  H3X 

t  A  *  I'  .  X  *  ■  #  i '  i  ♦  f  »  i  *  w .  'I*.  (  d  *  <  )  1  '  i  i 

_=0> 

(  *  •  •  •  *  *  •  f  ;  - 

IMX-  GC  •  !  )  m>3U*X  )*S03M  =  AU 
(OSH  MY-'  0  •  I  )  -  *S  \  *  Y)  *SHX  = 

soum  ox j+«  x  \ sHx*Y  +  r.  s  \  stoMx  )  »:>u 

030*00 T  +  SY.J*  SYS=aU 

(  /O  \  [  -  )  T  A=  (S  V  - 1  V  )  .  V  \  T  >H  a  (  0  v  - 1  V  )  T  >  =  (  '  V  -  f  V  )  s 

v  i  1 f  u  .  !  i  +■  {  vj  U  +  1  +■  3  +•■,,;)>.  ]  =  I  i  3 

(  fl0\G.  [-G  .1  )*M3X*9c  8?.  I  *fc 

IT*fl-HU3=IU  81 

IT  TA  0A0  OK.IHjQIciUCO  Y  i  IU  X03H3 
{03u+>iiAu>+J3i*3=en3au 

(  *iG  •.  i  i «  w_.  *  .»*  •  ♦  A •  t  C.  *  S'- J  «  J  X  *  I  I  )  V  U  /■  ;  t'  i  I  J  J  A 3 

(  Su0U*  (  X  —  OCU  C ■-  .1  )+C3  U*X  )  *SG3l'  =  AU 
(  u  a  .  Ml  -000^  0  .  I  )  +XH,  *  Y  )  *SHX=8  • 

XG;  M  Xj+,  G\SHX*Y+  ..  .  S\Sg0I-  *X  )  =3  ■ 

03‘J A  00 1*  +SYiU*SHHeQU 
3U  +  G/.U*  SOS+SuU*SOH=3. 
i *  (  GU+0U+OO+AG  )  +30  =  YI30G 

I T  =  I TM 

<0/  ,  llu,  j  ’  U  f  IT  )  H  .IT  J  J  A3 

(  IT\( IT-I TH)  ) 3 3 A  =  T  T 
VI  GT  00  { I  ,  .C.T0.TT  >31 
(G.  *TG.3i  ■  •  HO. S  •  T_  «Q> ) 3  I 

3?\ i  v=sv 

(  .AI*.AAi)  \H3X  *  *  T  *.GAei=IV 

G=0> 

Y^i^/i  T  *•  j  Ti^AlcuuG  I  h  ~  a!  h<j  i  0 ..  jXIhMC  3  'iTAdAlGA 

J  JG 71  30  YHUflTUB 
GS  OT  00  ( • A.TJ ,3X ) 3  I 
(II-  T ) *  I dO  •  C  +  {  I T \ S  T ) 0OJ A*f  ? •  '=0  33 

AS  OT  00 

'  • S**S T ) *. £45  > “ (1 T  .+  < i T \ST ) I 

(  •  G  G  C  I 

3 

S 1 \c  *  I =A  AS 
II  \  ^  •  I  =  d 

V- (  -  ■  *  .  ul  «  .  i  •  {  ,\J  •  -  )  *SODM*3»SO  : 

( ( * .S**I 

'  *  1  "  •  A  >  ♦•09S  ;  )  .  .  )  X  *  i  ... 

(  (  c  •  S  *  *  I 
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S02= (MOS+F* ( EXC+X*MCO/2.+Y*XH2/2. ))*(9.529*ALoG(T2/T1)+173. 0*2.0* 

1  (  A**Q. 5-B**Q. 5 ) -153  0. 0* ( A-B )  ) 

SH20  =  F*XH2*( l.-Y ) * ( 1 7 . 8  7*ALGG ( T2/T1 ) +5 9  7 . *2 . * ( A**0 . 5“B**0 • 5 ) “7 500 . 
1*( A-B) ) 

SN  2  =FN2  * ( 1 • 0  +  F ) *( 7 .48*ALUG (T2/Tl)+3470.*( A“B ) “5  8  0000. * ( A**2.~B**2 . 
1  )  ) 

SH2=F*XH2*Y* ( 3.77*ALUG( T2/T1 ) +0 . 0005 78* ( T2“T 1 ) “40 . * ( A**0 • 5“B**0 . 5 ) 
1  ) 

SUM=SFU+SC02+SC0+S02+SN2+SH20+SH2 
SU=  XCH  '•1.9859 *ALOG(CR  ) 

C 

PT2=T2 

26  CALL  TEMP  ( T2 *  SUM , SU , D I F * KO ) 

P2=CR*1 .0*T2/T1 

IF  (  KO.LT.  NOR.DIF.GT.O.O  )  GO  TO  20 
CALL  ENERGY  (  T2  > XC * UCQ2 , UCO , U02  ,  UN2 , UH20 , UH2  ,UFU  ) 

UA=MC02*UC02 
UB  =  XH2  *UH20 
UC=  EXC*U02 
UD=  PN2*UN2+MC  )*UCO 
UE=P02*UC2+PN2*UN2+UFU 
U2  =  F * (  JA+U8+UC+UD)+UE 
LS=U2+CEF 

CALL  C3USTN  ( XC » LS *  PC  * MC02 , XH2 , PN2 , EXC ,MCo , X3 , Y3 , T3 , RL , j2  ,  P2 , XCH , X 
1M  ) 

XC3=XRS+XM 

REVERS I BL  £  ADIABATIC  EXPANSION  -T3  TO  T4  -  AT  CONSTANT  ENTROPY 
X  AND  Y  BOTH  DECREASE  DUE  TO  TEMPERATURE  DROP 

ENTROPY  AT  T3 
PP  3  =  P3 

P3  =  XC3*T3*P2/ ( XCH*  T  2 ) 

27  CALL  ENTR0P(T3»X3»Y3*MC02*MC0,EXC*XH2*PN2,P3*XC3,S3) 

V3 =XC3*1 545. *T3 / (P3*14. 7*144. ) 

V4=CR*/3 

TO  FIND  T 4 

28  LS=0 • 0 
P  T4  =  P4 
KO  =  0 

L  L  =  1 8 

29  I F ( T4.GT. 3000. )  GO  TO  30 
X4  =  0  a  0 
Y4  =  0  •  0 
X  M  =  0  •  0 
T  T4=  T  4-1 0 • 

T4=TT4 
GO  TO  31 

30  KL=LL+1 

CALL  CBUSTN  ( XC » LS *  PC » MC02 > XH2 »XN2 » EXC *MCO , X4 * Y4 , T4 * NL , T1 ,P3»XCH*X 
1M) 

31  LL=KL 
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ST=STS 
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C  TO  FI ND  P  + 

XCH4=XRS+XM 
P4  =  XCH +  *  T  4 / ( XCH*T 1 ) 

PS4=S4 

CALL  ENTROP  (  T4,X4»  Y4  »MC02  »i^CO  ,  EXC  *  XH2  * P'^2  ,  ^4  , XCH4  , S4  ) 
I F ( S4.  GT • S3 )  GO  TO  29 
I T= 1 00 . * ( S3-S4 ) / (PS4-S4) 

K  T  =  T  4 

T4=  I  T  +  KT+ 1 0 

33  CALL  ENTR0P(T4*X*Y*MC02»MC0,EXC»XH2»PN2*P4*XCH4*S4) 
IF(T4.LT.1000. )  GO  TO  35 
KO  =  2 

CALL  TEMP(  T4*S4»S3»L»IF>KO) 

IF(DIF.GT.O.O)  GO  TO  33 
35  CONTINUE 

EXHAUST  PROCESS  AND  CLEARANCE  GAS 
V4  =  XCH4*1 5  45 •  *T4/ ( P4* 1 44 . * 1 4 . 7 ) 

DV= V4-V 1 

I F ( ABS(DV) .LT.50. )  GO  TO  32 

I SENTROP I C  EXPANSION  TO  T5  AT  14.7  PS  I 
ASSUMING  NO  DISSOCIATION  AT  T5 
32  P  5  =  1 .00 
PT5=T5 
X5  =  0  •  0 
Y  5  =  0  •  0 
T5=PT5+250. 

XCH5  =  XRS 
KO  =  0 

38  CALL  ENTR0P(T5»X5*Y5»MC02*MC0,EXC»XH2»PN2*P5»XCH5»S5) 
CALL  T  IMP ( T5  >  35  *  S3  *  D I F  *  KO ) 

IF(KO.LT.2.0R.DIF.GE.O.O)  GO  TO  38 
TT=ABS( (PT5-T5)/T5) 

I F ( TT.lT.O.OI )  GO  TO  39 
I  T  1  =  T 1 
T1=IT1+5Q 
GO  TO  17 

39  V5=XCH5*1545.*T5/ ( 14.7*144. ) 

FF  =  T  1  *  <CH/ (CR*T5*XCH5 ) 

DF= ABS ( ( F-FF ) /FF ) 

LF=  1 

GF  =  V2 / V5 
F  =  FF 

IFIDF.LT. 0.01)  GO  TO  40 
I  T  1  =  T  1 
T 1 = I T 1+50 
GO  TO  16 

TO  FIND  THERMAL  EFFIC.  AND  BMLP 

40  CALL  ENERGY  (T3*XC*UC02»UC0*UO2*UN2*UH2U*UH2*UFU) 
WRITE (6 *65)  PE*T1,T2»P2»T3*P3*T4,P4,X3*Y3,X4,Y4,F 
UA=MC02*( X3*UCO+( 1.0 0000- X3 ) *UC02 ) 

UB=XH2* ( Y3*UH2+( 1.00000~Y3 )*UH20) 

UC= ( X3  *MC02/2. 0  +  Y3*XH2/2.+EXC ) *U0  2 
UD=PN2*UN2+MC0*UC0 


f . 


'  I  • 


4  H  •  • 

x  +  ifl'.s^h  .• 


{  •  I  ...  A  *  1  •  '•  1  1  • 


l  K  i  '  )  \  t  l  4  1  J  X  = 

4  >3=aoh 

#4-  ♦  I  X  f  A  i  }  ' 1  O  At  >  J  jA 

. 

(  ^  )  \  (  A  -•  )*.  I  =  T  1 

AT  =  T> 
o  :  -  T >+  T  I  =  .  T 

*  *  »  *  Y  *  X  •  1  .  _ , 

:  ‘  0  ju  (  •  v.  •  .  ^  7  )  i  I 

s=ox 

(  o  •  -1 1  o  *  £  c  ♦  A  *  4'  T  )  H  •  i  T  J  J  A  0 
EE  Ci  Co  (  c  #  0  •  0  •  1  I Q  )  1 1 

iuHlTHGO  ?£ 


u  -j  jIA/ ■  y,  JO  un  A  t  j  JC  S  T  : U  Ai  1a  j 
..  #  .  A  A  I  C>M]  \  A  T  *  .eAc5I*AHDX  =  AV 


I  V-AV=V< 


( .  . t  j.  1 V'  ;  •  >1  . 

,  ; ; .  ■  h  1  ■  .  : . 

<:  T  T ,  Hu  i  ;  A I .  >c  c  I  <J  OH  OH  I MUSSA 

c..{=es  S£ 


e  T  =  c  T  H 
.  .  -  ■ '  X 

■- .  ^  =  e  ? 


.  ££+£TH=5T 
02X=eH0X 
C  =  0> 

It  •  -  HHf  i  *3X3fODM*  >Mt  <*  •  <  )  S( 


(  ■  0  '."HI  If  .  -  3  '  )  H.v:  T  J  J  A  ’ 

8 £  o T  Uo  (Of  .  f liU.MC. T J.OX ) 3 i 

(  ai  \  ( c?T  -c  Th  )  )  .  fj A  =  T  T 

(I  ..  J.JTHI 
I  T  =  I  T  I 


i  ■»!  T  I  =  I  T 
Vi  OT  00 


(  .  A-  I  -  .  :  )  V.  .  -  A  *ThD 

(  H0X*cM**  )\H0>*IT=11 
(  1  \ (  !  -  )  ) 88  A  =  3 

1=1  J 
?V\SV=10 
11=1 

OT  ..  !  Il.o.Tj.10)  II 

I  T  =  I  T  1 


oe  +  ITI=I  T 
di  OT  00 


0 

HjMc  UH  A  .01  J,  H  J  JAMRjHT  OH  11  OT  0 

■  f  •  ’  •  *  *  •••  •  »  ■  .'  •  i  .  J  A 

i«  f^Y*  AX#  *  ,  -  4  t  *.'*.*■  , 

(  'w(  k  (  1  :  j  0  .  D+.-O  • X  )  ^Svj0ri  =  A. 


(  o  .  1  '  (  -  .  '  Y  )  *  £  H  X  =  c 

0  *  (  ,  -  .  \  <  ■■  1  .  )  r. 


■  1  1  ,  ■ 


no  o  n  o 


87 


U3  =  U A+UB+UC+UD 

CALL  ENERGY  (T4>XC>UCU2>UC0«.U02>Ui'*2>UF)2o>uH2>0FU) 
U A  =  MCQ2  *  (  X4*UC0+  (1,00000- X4)  *UC02  ) 

UB=XH2* ( Y4*UH2+ ( 1 . G0000-Y4 ) *UFI20 ) 

UC  = ( X4*MC02/2.0 rY4*XH2/2 . +  EXC )*UJ2 

UD=PN2*UN2+ m c  0  * U  C  0 

U4=U A+UB+UC+UD 

W0RK=U3-U4+U1-U2 

ETA=100.*WORK/HEAT 

V 1  =  XCH*1 5  4  5 • *T  1  / ( 1  +.7*144.  ) 

V2=V1/CR 

BMEP  =  WORK*778./ ( 144.*( V1-V2  )  ) 

WRITE (6 >57)  U 1 >  U2  >  U3  >  U4  >  WORK  >  ET  A >  VI »  BMEP 

T  T  1  =  T 1 

TT2=T2 

TT4=T4 

TT  5  =  T  5 

DPC=PC+0. 10 

PC=DPC 

11  CONTINUE 
PCR=CR 
CR  =  PCR+ 1 • 0 

IF(CR.LE.ll.O)  30  TO  6 
GO  T 0 ( 2  >99 ) >  K I 

51  FORMAT (F6.0) 

52  FORMAT ( 1H1 >40X >8H  FUEL  IS>2X>A6) 

53  FORMAT ( 1HK*15X>13H  PERCENT  A1R=>F7.2) 

54  FORMAT ( 1HJ > 20X>20H  COMPRESSION  RATIO  =>F7.2) 

57  FORMAT ( 1 H  >8F16.2) 

6  0  FORMAT ( 1HJ>6F15.4»3F12«>4) 

65  FORMAT ( 1HK»9F9.2»5F10.6) 

99  STOP 
END 


SIBFTC  ENTROP  DECK 

SUBROUTINE  ENTROP ( T 3 > X > Y  > MC02 , MCO , EXC > XH2 ,HN2*P,XC*S) 

REAL  MC02  >  MCO 
T  =  5  2  0  *  0 
A= 1 • 0/ T  3 
B=1 •0/520*0 

I F ( T3.LE.600.  )  WKlTE(6*59)  T 3 , X > Y > MC02 ,MCo , EXC , XH2 > P > XC 
T  LuG  =  A  LOG ( T3/T  ) 

SCo2=Mc02* ( 1 ,-X)*( 16.2*TLug+6530.*( A“B ) “705000 .* ( A**2 . 0“B**2 . 0 ) ) 
SCO= (MC0+MC02*X)*( 9.46*TLOG+3290.*( A“B )“535000.*( A**2.0“B**2.0) ) 

SO 2= (EXC+X*MC02/2.  +  Y*XH2/2. ) * ( 1 1 . 5 1 5*TLOG+ 1 73 . 0*2 . 0* ( A**0 . 5 “3**0  •  5 
1 ) - 1 5  3  3 • * ( A-B )  ) 

SH20=XH2*( l.-Y) *( 19. 8  6*TLOG  +  597. *2.* ( A**0. 5“B**0. 5 ) “7  5 00 . * ( A“B ) ) 
SN2  =  PN2*(  9.47*TLOg  (-3470  .*(  A“B  )  -58  000  0  .*(  A**2  .  “B**2  .  )  ) 

SH2  =  XH2*Y* ( 5 .76*TLOG+0.  0005 78* ( T 3“T ) ~40 . * ( A**0 . 5“B**0 . 5 )  ) 

I F ( T3.LE.4000. )  GO  TO  27 

EH2=SH2~0. 00033* ( ( T3“T) ~4000.*TLOG)*XH2*Y 
SH2=EH2 

I  F  (  T3.LE.5000.  )  GO  TO  'J 

EO2  =  SO2  +  0.0CC0  5*(  (  T3“T ) “4000 . +TLOG ) * (  E  XC  +  X*i/iC02  /2  .  +  Y*XH2  /  2  .  ) 
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S02=E02 

27  SUM=SC02+SC0+S02+SN2+SH20+SH2 

TO  FIND  PARTIAL  PRESSURES  OF  CONSTITUENTS  AT  T 3 
A= (  1 ,0-X) *MC02 
B= ( EXC+X*MC02/2.+Y*XH2/2. ) 

C= (  1 .00-Y ) *XH2 
D=  (  MCO  t-MC02*X  ) 

E  =  PN2 
G=  XH2*Y 
AP=A*P/XC 
BP  =  B*P / XC 
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C 


21 

22 


23 

24 


25 

26 


28 

29 


CP=C*P/XC 
DP=D*P / XC 
EP=E*P/XC 
GP=G*P/XC 

SUM  OF  1  •  98  5  9 ( M ) LN ( PP ) 

IF ( AP.LE. 0.0 )  GO  TO  21 
SA=A*ALOG( AP ) 

GO  TO  22 
SA  =  0 • 0 

IF ( CP.LE. 0.0 )  GO  TO  23 
SC=  C*A  -OG ( CP ) 

GO  TO  24 
SC=0.0 

I  F ( EP. lE. 0 .0 )  GO  TO  25 
SE=E*ALOG ( EP ) 

GO  TO  26 
SE  =  0 . 0 

IF(BP.LE.O.O)  GO  TO  28 
SB=B*A -OG ( BP ) 

GO  TO  29 
SB=  0 • 0 

I F ( DP.LE.0.0)  GO  TO  31 
SD=  D*ALOG ( DP ) 

GO  TO  32 


31  SD=C • 0 

32  I F ( GP.LE.0.0)  GO  TO  33 
SG=G*ALOG( GP ) 

GO  TO  34 

33  SG  =  0 • 0 

34  SUS=1.9859*( SA+SB+SC+SD+SE+SG) 

SD I S= ( MC02*X*121 181 .  +  XH2*Y*10348  6.  ) / T  3 


S  =  SUM-SUS  +  SD I S 
RETURN 

59  FORMAT  UH  *8F16.2) 
END 
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$  I  BFTC  ENERGY  DECK 

SUBROL  riNE  ENERGY ( T  >  XC  >  UCG2 > UCO , U02 , UN2 , UH20 , UH2 , UFU ) 

C  UT-U520=INTEGRAL  OF  CV*DT  =  INT ( CV* ( T~520  )  )  CV=CP~1.9859 

TE=T-520. 0 

r= i  •  o / r 

S=1 . 0/520. 0 
V=  T  **2 • 0 
W=  5 2  0 • J**2.0 
X=  T**0 • 5 
Y= 52  0 • 1**0 • 5 

UC02=( 16.2-1.°859)*TE_6530.*ALOG( T/520. ) -141.* ( K-S ) *10000. 

UCO= ( 9 .46-1.98  59 )*TE“32 90. *AlOG( T/520-  )~107.*(R-S)*10000. 

U02= ( 11.515-1.9859)*TE-173.*2.*(X-Y ) +1 5 30 . *ALOG ( T / 520 . ) 

UN2= ( 9.47-1.9859 ) *TE“3470.*ALjG( T/520. )-116.*(R“S ) *10000. 

UH2  0= (19.86-1.9859 ) *TE~597 . *2 .* ( X~Y ) +7500.*ALOG (T/520. ) 

UH2  = ( 5.76-1.98  59 ) *TE  +  0. 578* ( V-W ) / 2 00  0. +  40 .* ( X~Y ) 

I F ( XC.LT.4.0)  GO  TO  1 

UFU= ( 7.92-1 .98  59 ) *TE+0. 0601* ( V-W ) / 2 . 0 
GO  TO  3 

1  UFU* ( 2.258-1.9859) *TE+. 032* ( V-W) /2. -543.* ( ( T/100. ) **3 .~5 .2**3 . ) /30 
10. 

3  IF ( T.LE.4000. )  GO  TO  2 

UH=UH2~.33*( . 5*T**2.-4000.*T ) /1000. 

UH2=UH 

IF ( T.LE.5000. )  GO  TO  2 
UO=UO2+.05*( .5*T**2.-4000.*T ) /1000. 

U02  =  UO 
2  RETURN 
END 


S I BFTC  TEMP  )ECK 

SUBROUTINE  TEMP ( T , SUM , $U , D I F » KO ) 

D I  F  =  SUM-S'J 
TE  =  T 

IF(KO.GT.O)  GO  TO  2 
I F ( DIF. LT. 0.0)  GO  TO  2 
I T  = ( T  +  5.C ) /10. 

TE  =  IT*  10 
DT  =  -100 • 

GO  TO  3 

2  I F ( KO. iQ. 1 • AND. D IF.GT .0.0)  GO  TO  5 
I F ( KO. EQ. 2 )  GO  TO  6 
KO  =  1 
DT  =  1 0 • 

D I  =  D  I  F 
GO  TO  3 

5  KO  =  2 

6  DT  =  - 1 • DO 
D  I  =  D  I  F 

8  T  =  T  E  +  DT 
RETURN 
END 
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$  I  BFTC  CBUSTN  NODECK 

SUBROUTINE  CBUSTN  ( XC,LS  »PC  ,MCu2 ,XH2 ,HN2 , EXC ,MCo , XX > Y Y , T  * Nl , T2 , P2 , 
1 XCH  >  XM ) 

C 

C  DISSOCIATION  OF  002  AND  H20  ONLY 
C 

DIMENSION  TEPR ( 2  5 ) *  02 ( 25 ) * N2 ( 25 ) >C02(25)  ,H20(25)  » H2 ( 50 )  »C0 ( 25 )  > 

1 LKC02 ( 2  5 )  »  LKWG (25)  »KC02(25)  ,KTC02(25)  >X(25),Y(25) 

REAL  N2*LKC02,LKWG»MR,KCO2»KTC02*KWG»LS»M02 *MC0*MC02 *MOS 

INTEGER  TEPR 

I F ( LS.  IQ. 0.0 )  GO  TO  4 

IF(KL.GT.O)  GO  TO  2 

DO  3  1*1*23 

READ (5 >20)  TEPR( I ) >LKC02( I ) >LKWG( I ) 

KL=  1 

3  CONTINUE 
2  I  L  =  1 

1  1=  I  L 

C  CHEMICAL  ENERGY  OF  DISSOC  OF  H20  AND  C02 
CEH2=103486.*XH2 
CEC0=121181.*MC02 
C 

C  SUMMATION  OF  TERMS  IN  ENERGY  EQN • 

TEMPR=TEPR ( I ) 

CALL  ENERGY  ( TEMPR,XC>C02 ( I ) *CO(I)  *02 <  I  >  >N2 ( I ) *H20(I)  ,H2 ( I ) *UFU) 

IF ( EXC.GT. 0.001 )  GO  TO  30 
U02  =  0 . 000  0 
GO  TO  31 

30  U02=EXC*02 ( I ) 

31  IF(MCO.GT.O.OOl)  GO  TO  32 
UCO=0. 000000 

GO  TO  33 

32  UCO  =  MCO*CO (  I  ) 

33  SUMN  =  MC02*C02 (  I ) +XH2*H20 (  I ) +PN2*N2 (  I  ) +EXC*02 (  I ) +U02+UCO 
SUMX=MC02*C0( I )-MC02*C02( I ) + ( MC02 /2 . 0 ) *02 ( I ) +CECO 

SUM Y=XH2*H2 ( I )-XH2*H20{ I )+( XH2/2. ) *02 ( I )+CEH2 
A=LS~SUMN 

I F ( A.GE.O. 001 )  GO  TO  34 
I L= I +1 
GO  TO  1 

34  B  =  A / SUM  Y 

35  C=SUMX/SUM Y 

c  obtain  log  kvg  from  tables 

K WG= 1 0 , OCOO**LKWG ( I ) 

D=KWG-1 .000 

C  KWG=X( 1 - Y ) / Y ( 1~X)  SO  Y=X/ ( KWG~ ( KWG"1 ) X )  SO  Y=X/ ( KWG~D*X ) 

C 

C  EQUATE  Y@S  (  B-C-v-X)*(  KWG-D*X)  =X 
CC  =  B*K WG 
AA  =  C*D 

BB  =  B*D+C*K WG+ 1 • 000 

DD=BB* ^2-4.000*AA*CC 

X( I )=(BB-SQRT(DD) )/<2.000*AA) 

Y (  I ) =  X  C I ) / ( KWG~D*X ( I  )  ) 

C  TO  FIND  KC02=  (  (  P(?C02  )  **2  )  *  (  P@02  )  /  (  P@C02  )  **2 
C  NO  OF  MOLES  OF  PRODUCTS  =MP 
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MP  =  MC02  +  XH2  +  MC0+EXC  +  IjN2+  (MC02*XU)+XH2*Y(I)  )  /2. 0 
R  =  X (  I ) **2 

S= ( 1 .0000-X( I ) )**2 

W=(MC0i/2.  )  *X  {  I  )+(XH2/2*0)*Y(  I  ) 

TEMPER- TEPR (  I  ) 

KCQ2 ( I ) =( R* TEMPER* W*P2 ) / ( S*XCH*T2 ) 

KTC02 (  I  )  =  lu.oG0**LKC02 ( I  ) 

I L= I +1 

I  F (  I • EQ • 2  3  )  GO  TO  12 
I F ( LS.GT.0.0 )  GO  TO  6 
4  I  =  KL 

KO  =  0 
XX=0 • 06 
WX=7000.*XX 
KWG= 1 0 • Q0**LKWG ( I ) 

D=KWG-1 .00 

KTC02( I )=10.00**LKC02 ( I ) 

T  =  TE  PR (  I ) 

5  YY=XX/ (KWG-D*:;x  ) 

IF( XX. 3T.0.1 )  GO  TO  14 
IF ( XX. GE. 0.00015 )  GO  TO  11 
XX=0.0001 
YY=0.C3001 
11  A=XX**2.0 

B= ( MC02*XX+XH7*YY ) /2 . 0 

D= (  l.C-XX)**2.0 

KC02 ( I ) =B*  T*A / (D*T2*XCH) 

CK=  1 00 30000 , *KC02 ( I ) 

TK=1000000  0,*r.TC02  (  I  ) 

CALL  TEMPI WX»CK,TK,DIF»KO) 

XX=WX/ 7000 • 

I F ( KO.LT.2.0R.DIF.GT.O.O )  GO  TO  5 
GO  TO  14 

6  I F ( KC02 (I) .GT. KTC02 ( I ) )  GO  TO  8 
IF ( T.LE.3000. )  GO  TO  14 
P  T  =  T  EPR ( I  ) 

PX=  X (  I  ) 

P  Y=  Y  (  I  ) 

PK=  KC02 ( I  ) 

PKT=KTC02 ( I ) 

P  K  w  =  K  W  G 
PMP=MP 
GO  TO  1 

8  F  =  KC02 (  I ) -KTC02 ( I ) 

WR I T  E ( 6  » 2 1 )  TEPR( I  ) *SUMN,SUMx*SUMY»MH,KCu2(  1 )  »KTC^2(  I  )  ,X( I )  >Y(  I  ) 
IF ( I • GT • 1 )  GO  TO  10 
K  L  =  2 
T=5400. 

XX=X ( I ) 

Y Y=  Y  (  I  ) 

WR I TE ( 6  »27 ) 

GO  TO  14 
10  G=  PK  T-PK 
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IF ( TEP  *( I ) . LT.  5000 )  GO  TO  7 
SPA=200. 

GO  TO  9 
7  SPA=100. 

9  I T=  F*S  PA/ ( F  +  G ) 

DT=  IT 
T  T=  TEPR ( I  ) 

T  =  T  T  +  DT  +  5 . 0 
U=  X (  I ) -PX 

XX=PX+( SPA-DT )*U/SPA 
V=  Y (  I ) -PY 

YY=PY+( SPA-DT ) *V/SP A 

dmp=pmp-mp 

XMP=MP+DMP*DT/SPA 
MP  =  XMP 
GO  TO  14 
12  WR I TE ( 6  »  2  5  ) 

KL  =  2 

14  XM= (MCD2*XX+XH2*YY ) / 2  •  0 
RETURN 

20  FORMAT ( I4»F6.2*F7.3) 

21  FORMAT ( 1HJ* I5,3F12.0*F8.3,4F14.7) 

25  FORMAT  (  1HJ  »44H  TEMPERATURE  BELOW  3000(?>R  MORE  DATA  REQUIRED) 
27  FORMAT ( 1HJ » 14H  T3  ABOVE  5400) 

END 
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APPENDIX  D 

SAMPLE  CALCULATION  OF  TEST  RESULTS 


Test  No,  5,  Run  43 

C .  R.  =  9,77 

Spark  20° 

Revs  2333 
Time  3 . 845 
Fuel  Temp,  7  8°F 

Fuel  Pressure  2  x  0.68  =  1.32"H20 
Air  Temp,  83 °F 

Air-flow  Manometer  0.227  in.I^O 


Barometric  Pressure  =  27,590 

in.Hg. 

Cylinder  Temp.  585  °F 
%  C02  8.4 
%  02  7.4 

%  0  0,0 

Dyna.  Load  10.9  lbs. 

Fric.  Load  10.1  lbs. 


Gas  Pressure? 


1,32  (29.92) 

33.899  (12) 


0.0971  in, Hg . 


Saturation  Pressure 


P 


.3631  +  8(o 01438) 


.3631  +  .1150 

(0 . 4781) 

27.590  +  0.971  - 


29.92 

147696 


29.  92 
14.696 


0.973  in.Hg. 


0.973  =  26.714  in  Hg 


0.500  cu.  ft.  of  gas  in  3.845  mins. 
Correct  for  Barometric  conditions. 
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Vol . 


. 5  (26 . 7 14) 520 

29. 92 (460+78) 


431  ft" 


Wt, 


.431(60) 

8.503 (3.845) 


=  0.791  lb/hr. 


Correction  factor  (CF) 


r.p.m, 


Torque  = 


29 

.90 

27 

.590 

29 

.90 

27 

.59 

2333 

3. 

845 

10 

.52 

460  +  78 


545 


0.989 


1.078 


606  r.p.m, 


12 


(10.9) 


9.56  ft-lb. 


BHP 


2?r(9. 56)606 

33,000 


1.103  HP 


FHP 


2t r 


10 . 52  \  10.1(606) 


12 


33,000 


1.022  HP 


CIHP  =  1.078  (1.103  +  1.022) 

=  1.078  (2.125)  =  2.29  HP 
C.B.H.P.  =  2.29  -  1.022  =  1.27  K 


Thermal  Effic.(r)) 


2.29(2545)100 

19929 ( . 791) 


37.9% 


Brake  Thermal  Effic.  = 


1.27 (2545)100 

19929 (.791) 


20.5% 
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Meclu  Efficiency 


Fuel/ CIHP 


Fuel/BHP 


Air  Density 


Head  of  air 


Velocity 


Flow  of  air 


Fuel-air  Ratio 


~~  (100)  =  55.5% 

~y|  =  0.346  lb/ CIHP 


.791 

1.27 


0.623  lb/BHP 


27.590(14.696)144  3 

29 1 92  (5381  53.34  =  °0679  lb/ft 


.227  (62.4) 

.0679  (12) 


17.4  ft.  of  air 


0.62  y  64.4(17.4) 

0.62(33.45)  =  20.75  f.p.s. 


20.75  ir  ( .  9983 )  2 

4(144) 


60  ( o  067  9) 


1.078 


20.75  ir(.9966)60(.0679) 

4(144) 


1.078 


0.495  lb/min. 


0.791 

0.495 (60) 


0.0266  lb. fuel/lb. air 


1 


Air-fuel  Ratio 


.0266 


37.6  lb. air/lb. fuel 
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B„m„  e. p„ 


I oitio  e„  p. 


33,000(1.27)4(2) 

ir(3. 25)  2  („375)606 


33,000(2.29)4(2) 

tt(3.25)2  (.375)606 


80.1  psi 
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APPENDIX  E 

PROGRAM  for  CALCULATION  of  TEST  DATA 


$UOB 

525075 

G.R.  POND 

$T  IME 

5,3000 

$  I BUOB 

POND 

$  I BFTC 

POND 

DECK 

C  CALCULATIONS  FOR  C.F.R.  ENGINE  TEST 

DIMENSION  Run (50)  *REVS( 50 ) *TIME< 5 0 ) , temp (50 ) *DP(50) ,TEM(50) *Co2(50 
1 )  *02 ( 50 ) ,CO( 50 ) *  LOAD (50)*FRIC(50)  *RPM (50)»LB(50)*  ToR ( 5  0 )  *BHP ( 5  0 )  »E 
2TA ( 50  )  *  FHP ( 50 ) » IHP ( 50 ) *MEFF( 50) *FU( 50 ) »RATIU( 50  )  *Z ( 50 ) »E( 50 ) *FA( 50 
3 )  »  BMEP ( 50 )  > A IR ( 50 ) *CBHP ( 50  )  »CMEP (50)  *FUEL(20)»SET(50)  ,CF(50>»CK(2 
4u)  *  BAPRES  (  20)  ,  TEST  (20)  ,SPARM50)  ,mmucM50)  *  FUE  (  20  )  >  TE  (  5  0  )  » AF  (  5  0  )  » 
5HYC(  50)  >FUA(  50)  >COM(  5  0  )  >Men  ( 50  >  * TS ( 50 )  , FU I  HP ( 50 ) * CYTEM ( 50  >  *FUTEM  (  5 
60 )  *  PRE  3 ( 5 0 ) 

REAL  LOAD,  LB  *  I  HP , N2 » MASS » MEF F »M W , MEN 
INTEGER  RUN, NO, I 
DD  =  5  00  •  0 

READ (5,99)  M,  PROPA 
DO  24  J  =  1 »  M 

READ(5,100)BAPRES(U) , CR ( J ) , TEST ( J ) ,NO,FRlC( J) ,FUEL(U) ,FUE( J) 

WRITE (6, 104) 

WRITE(6»105)TEST(J) ,CR(U) ,FUEL( J) ,FUE(d) ,BApRES(J) ,NO 
DO  1  1=1, NO 

READ(5»106)RUN(I),SPARK(I),REVS(I),  TIME  (  1  )  ,  FUTEM(  1 ) » PRES (  I  )  ,  TEMP  (  I 
1 )  , DP (  I )  , TEM(  I )  ,CU2 ( I )  ,02 l I )  ,C^( I ) , LuAD(  I  ) ,CuM( I  ) ,MEn(  I )  ,  TS (  I  ) 

RPM (  I ) =REV S (  I  ) / T I ME (  I  ) 

1  CONTINJE 

I F ( FUEL ( J ) • EO,  PROPA)  Gu  TO  3 
WR I TE ( 6 , 1 02 ) 

DO  2  1=1, NO 

2  WRITE(6»111)  RUN ( I ) , SPARK (  I  )  ,REVS( I )  , T I M E ( 1 ) , TEMP(  I ) , DP ( I ) , TEM ( I )  , 
1C02 ( I ) ,02 ( I ) »CO( I ) , LOAD ( I ) ,COM( I ) , MEN ( I ) ,TS( I ) 

GO  TO  5 

3  WRITE(6»1Q3) 

DO  4  1=1, NO 

4  WRITE (6, 112)  RUN ( I ) , SPARK (  I  )  »REVS( 1 )  , T I ME ( I ) , FUTEM ( I )  ,  pReS (  I ) , TEMP 
1 ( I  )  *DP ( I ) , TEM( I ) ,C02 ( I ) ,02 (  I ) ,CO( I ) , L^AD ( I )  ,CdM (  I  ) »MEn (  I  )  ,  TS (  I  ) 

5  CONTINUE 

WR I T  E (  3,113)  FRIC(U) 

PR  ESS= ( BAPRES (U) *14.6  96 ) /29.92 
IF ( FUEL(U) .EO,  PROPA)  GO  TO  7 
C  FUEL  IS  GASOLINE 
X=  7  • 

Y=  1 6  • 

H V  =  1 9 1 37. 

MW  =  100. 198 
ST0IC=  3.068 
DO  6  1=1, NO 

C  FUEL  CONSJMPTION  (50  GRAMS) 

LB ( I )=( 0.050*2.2046*60. ) /TIME ( I ) 

6  CONTINUE 
GO  TO  11 
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C  FUEL  IS  PROPANE  GAS 

7  X=  3  • 

Y=  8  • 

HV  =  19929* 

MW=44 • 098 
STO I C=  0 • 064 
DO  10  1=1, NO 

C  FUEL  CONSUMPTION  IS  1/2  CUBIC  FOOT 

GAPR=2.*PRES(I) *2 9. 92/(33. 899*12.0) 

I F ( FUT £M( I ) .GT.70. )  GO  TO  8 

PS AT= ( 0.3631“ ( 70. -FUTEM( I))*0. 010 68 >*29.92/14. 6 96 
GO  TO  9 

8  PSAT=( 0.3631+ ( FUT EM ( I ) -7 0 . ) *0 . 0 1 438 ) *2 9 . 92 / 1 4 . 6 96 

9  PT=BAPRES( J) +GAPR-PSAT 

VOL  =  0.5*PT*52f)./  (  29 . 92*  (  FUTEM  (  I  )+459.9)  ) 

LB (  I  )=30.0*VOL/(8.503*TIME( I )  ) 

10  CONTINUE 

11  CONTINUE 

DO  13  1=1, NO 

C  CORRECT  ION  FACTOR  (TO  29.90  IN  FIG.  AND  8  5DEGREES  F.) 

CF (  I  )  =  ( 29.90 /BAPRES(J)  )*SuRT(  ( 45  9. 9+ TEMP ( I )  )/ 545.0) 

C  TORQUE 

TOR (  I  )  =10. 52*LOAD(  I  )  / 1 2 . 0 
C  BRAKE  HP 

BH P (  I  )=(2.0*3.14159*TOR( I  )*RPM( I )  )/ 33000.0 
C  FRICTION  HP 

FHP(  I  )  =  (2.*3.14159*RPM( I)*FRIC(U)*10.52)/(33000.*12.) 
C  INDICATED  HP 

I  HP (  I  )  =  ( BHP (  I  ) +FHP ( I  )  )*CF(  I  ) 

C  CORRECTED  BHP 

CBHP ( I ) = I  HP ( I ) -FHP (  I  ) 

C  INDICATED  THERMAL  EFFICIENCY 

ETA (  I  )  * ( I  HP (  I  >*2  545.0*10  0.0)/ (HV*LB(  I )  ) 

C  MECHAN I CA_  EFFICIENCY 

ME  FF (  I  )  =  ( CBHP ( I )*100. )/I HP ( I ) 

C  FUEL  PER  I  HP 

FU I  HP ( I ) =LB ( I ) / I  HP ( I  ) 

C  FUEL  CONSUMPTION  PER  CBHP 
FU  (  I  )  =  L B (  I  )  /CBHP (  I  ) 

C 

C  **FUEL  TO  AIR  RATIO** 

C  MEASURED  (  CORRECTION  FACTOR  APPLIED  ) 

DENS=( PRESS* 144.0) / ( 53.34* ( TEMP ( I ) +459. 7 ) ) 

HE AD= ( DP ( I ) *62 . 4 ) / (DENS*12.0) 

VEL=u.62*( SQRT (2.*32.2*HEAD) ) 

FLOW= ( VEL* 3. 14259* ( 0 .9983**2 ) *1 5 . *DEN S/ 1 44 . )*CF( I ) 

RAT  1 0 (  I  )=LB( I )/(FLOW*60.0) 

C  AIR  TO  FUEL  R\TlO 

AF (  I  )=  -LOW*60./LB(  I ) 

C  FROM  ORSAT  ANALYSIS 
H2=0.51*CO( I ) 

CH4=0. 22 

TOT AL  =  C02 (  I  )  +02  (  I  )+CO ( I )+H2+CH4 

N2=lOO.OOO-TOTAL 

OX Y  =  N2  '3.76 
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C  ZCXHY+0XY(02) +N2 (N2 )  =  C02 ( C02 ) +  C0 ) +  0 . 3 ( CH4 ) +02 ( 02 ) +H2 ( H2 ) +N2 ( M2 ) +  2H20 
C  CARBON, OXYGEN, BALANCE 

Z(  I  )  =  ( C02 (  I )+C0( I ) +0 • 22 ) / X 
MASS= (0XY+N2 ) *28 .967 
FA (  I  )  =  ( Z (  I ) *MW ) /MASS 
e ( i )  =  ( fa ( i ) -Rat i o ( i ) ) *100.0 /Rat io ( i ) 

C  TO  CALCULATE  HYDROGEN  CARBON  RATIO 
XX=X*Z ( I ) 

W  =  2 .00* (0XY-C02 ( I ) — C 0(1) /2.0“0.22 ) 

Y Y  =  0  •  8  3+H2*2 .0  +  2.0* W 
HYC (  I  )  = YY/ ( 12. 0*X X ) 

FUA (  I  )  =  ( 12 .0*XX+YY ) /MASS 
C 

C  TO  CALCULATE  BRAKE  MEAN  EFFECTIVE  PRESSURE 

A  L  N  C  = (  (3. 14159*3.25**2 ) *0.37 5*RPM( 1  )  ) / 8 . 0 
BMEP(  I  )  =  ( 3  3000. 0*BHP (  I )  )/ALNC 
CMEP(I)=3300G.*IHP(I ) / ALNC 

13  CONTINUE 
C 

C  **  TO  SuRT  INTO  DESCENDING  VALUES  OF  FUEL  AIR  RATIO** 

L I M  =  NO- 1 

14  KNT  =  1 

C  INITIALIZE  KNT  IH  EVENT  ALL  RATIOS  IN  ORuER 
DO  15  1=1, LIM 

I F ( RAT  I  0 (  1  +  1 } .LT.RAT I  0 ( I )  )  GO  TO  15 
C  IF  RAT  1 0 (  1  +  1 ) .GT. RATIO!  I  )  INTERCHANGE  THEM 
S0RT=RATI 0 ( 1+7 ) 

S A  =  RUN {  1  +  1 ) 

SB=RPM { 1+1 ) 

SC=BHP ( 1+1 ) 

SD  =  FHP (  1  +  1 ) 

SE  = I  HP (  1  +  1 ) 

SF  =  TOR (  1  +  1 ) 

SG=  E  T  A (  1  +  1 ) 

SH  =  MEFF ( 1  +  1 ) 

S I  =  T  EM (  1  +  1  ) 

S J  =  FU (  1+1 ) 

SK  =  FU I  HP (  1+1 ) 

SL  =  CMEP( 1  +  1  ) 

SM  =  SPARK(  1  +  1  ) 

SN=CF ( 1+1 ) 

S0  =  CBHP ( I +1 ) 

SP  =  LOAD ( I +1  ) 

SQ  =  COM (  1  +  1  ) 

SQA  =  ME  4(1+1) 

SQB=  TS (  1  +  1  ) 

SR  =  HYC (  1  +  1  ) 

SS  =  FUA (  1  +  1  ) 

ST  =  F  A ( 1+1 ) 

SU=Z ( 1+1 ) 

SV=E ( 1+1 ) 

SW= AF ( 1+1 ) 

SX  =  C02 (  1  +  1  ) 
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*  •  J  /•  \  (  (  *  )  HH8*v.C  ?  f  )  =  (  I  )  HJM 

H  JA\(  1  UHI*.  vf>  =  (  I  IHjMD 

3UM I T  W03  £X 

#A'  i  im/i  Mi-,  -j  j  U  3  U  c  j  o  w ,  V  u » ,  i  u  *  i  J  36  i  t  vW'I  T  p  u  3  uT  *  * 

X  ~OH  =  M I J 

X  =  TH>  AI 

r:  j  o  /  u  ,.i  .u  i  T  A/  -j  j  A  ::i  Th^i  3SIJAITI/.I 

mij* x=i  ex  oo 
dl  >T  Oj  (  (  i  )  u  1  T  A  P  .  T  J  .  (  X  +  I  )  0  1  T  A  fl  )  3  I 
Jt  T  -OYIA-OM JTn i  (  i  UITAm.TO*  (  f  +  i  )01  TAfl  31 

( ;  +  i )c itap=tpo3 

( X  +  I ; MUP=A3 
(  X  +  I  )hHP=83 
( x+i ; HHP=oe 
(  X  +  I  )SH3=G3 
(  X  +  I  )  PHI =33 
( X  +  I  ) POT  =  33 
(X  +  I  I  AT  3  =03 
( X  +  1 ) 333M  *H8 
( X  +  I  IM3T  =  13 
( X+I  ) U  3  =  0  3 
(  X  +  I  )MHIU3  =  >3 
( X  +  I ) M  3  M  3  =  J  3 
(  X+  I  MflAP3*M3 
(  X+I  )  33  =  /s3 
(  X  +  I  )  H  H  8  3  =  0  3 
(X  +  i  HJA0J  =  M3 
(  i  +  I  ) M 03  =  03 
( X+I H  3M*A03 
( X  +  I  ) 3T  =803 
( X+I ) 3YH=P3 
( X+I ) AJ3=33 
( X+i  )  A3  =  T3 
4  (1+1)1*03 

( X+I ) 3=V3 
( X+I  )  3  A  =  W  3 
( X+I ) S03*X3 


3 


3 


3 


3 

3 
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S Y  =  02 (  1+1 ) 

SZ  =  C0 (  1+1 ) 

SZA=LB ( I ) 

RAT  10 (  1+1 ) =RAT 1 0 (  I ) 
RUN (  1  +  1  ) =RUN (  I  ) 

RPM (  1  +  1 ) =  R  PM ( I  ) 

BHP (  1  +  1 ) =BHP (  I  ) 

FH  P (  1  +  1 )  =  F  HP ( I  ) 

I H  P (  1  +  1 )  =  I  HP (  I  ) 

TOR ( 1+1 ) =TOR ( I ) 

ETA (  1  +  1 ) =E  TA (  I  ) 

MEFF (  1  +  1 ) =MEFF ( I ) 
TEM (  1  +  1 )  =  T EM ( I  ) 

FU (  1  +  1 ) =FU (  I  ) 

FU I  HP (  1+1 ) =  FU I  HP (  I  ) 
CMEP (  1  +  1 ) =CMEP (  I ) 
SPARK( 1+1 ) =SP ARK ( I ) 
CF (  1  +  1 ) =CF (  I ) 

CBHP ( I fl ) =  CBHP (  I ) 
LOAD ( 1+1 ) =LOAD( I ) 
COM (  1  +  1 ) =COM ( I  ) 

MEN (  1  +  1 ) =MEN (  I  ) 

TS (  1  +  1 ) =T  S (  I ) 

HYC (  1  +  1 )=HYC (  I  ) 

FUA (  1  +  1  )  =  FUA ( I  ) 

FA (  1  +  1 ) =F A (  I  ) 

Z(  1+1 )=Z(  f  ) 

E( I +1 ) =E ( I ) 

AF (  1  +  1 ) =AF (  I  ) 

C02 (  1  +  1 ) =C02 ( I  ) 
02(1+1) =02 ( I ) 

C0( 1+1 ) =C0 ( I ) 

LB (  1  +  1 ) =L  3 l I ) 

RAT  1 0 (  I  )  =SORT 
RUN (  I  ) =  SA 
RPM (  I  ) =  SB 
BHP (  I  ) ^SC 
FHP (  I ) =  SD 
I  HP (  I ) =SE 
TOR (  I  )  =SF 
ETA ( I ) =SG 
MEFF ( I ) =SH 
TEM (  I ) =  S I 
FU (  I  )  =  5 J 
FU  IHP ( I  ) =SK 
CMEP (  I ) =S  L 
SPARK ( I ) =SM 
CF (  I  )  =  5N 
CBHP ( I ) =S0 
LOAD ( I ) =SP 
COM (  I ) =  SO 
MEN (  I  )  =SQ A 
TS ( I ) =SQB 
HYC (  I  ) =SR 


(  f  + 1  )  =  ■- : 

( r  + 1)  dxu 
i  i  ;  j  a  s  . 

(  i  J  0  A T  A  P  = (I  +  I  )  u  I  TAP 
(  I  )  A  •  F: -  (  :  +  I  )  HUP 
(  )  --  (  +  1  )  h  HP 

(  UH'  =  (  i  +  I  )  SH  \ 
(  I  )  Hh  1  =  (  j  h- I  )SHi 
(  I  )  Hi;  i  =  (  i  +  I  )  SH  1 
(  I  )floT  =  <  I  +  I  )>UT 
(  I  )  A  T  j  =  (  I  +  l  )  A  T  i 
(  I  )  3  i  =  (  [  +  i  )  3  3  3f- 
(  I  ) M J  T  =  (  [  +  1  ) M3T 
(  I  )  L  X  I  +  I  )  J  3 
<  1  )S'K03=(  I+I  )HHIi  1 
(  1  )  H  j :  '  =  (  r  +  i  )  H  .  O 
'  is  .-''  (  Hi  )  •/- 

(I  )  3  ,  -  i  I  +  I  )  3  0 

(  i  )HH  ,J=(  ;-4  i  )HhJ) 

(  I  );;mUJ=[  I  +  I  )  0 AU J 
(  I  )f-XO*<  I+I  )  ;0D 
(  I  )  /i  :=(;  +  !  )  h  j t, 

(  I  )  A  T  =  (  I  + 1  )  £  r 
(I )  D  I*  1 1 +1  )  D  YH 
(  I  ) AU  3  =  (  I  +  I  ) AU3 
( I ) A3  = ( r+I ) A3 
(  ]  >X  I+I  )  + 

(  I  )  3= (  1+  I  ) 
(D  -  -1I+IHA 
(  I  )  SOD= ( I  +  I  ) SOD 
(  I  )  5sI0*  <  1  +  I  )  SC 
(IK  >  (  I  +  I  )  UD 
( I ) r J= ( I  +  I  >  3J 
0a*( 3  JOI TAP 
A  1= (  I  } i  UP 
X  (  I  )  MSP 
"••  =•(  I  )  SH  r 
0= (  I  ) S  H  3 
■Us  (  I  )  SHI 
H  (  I  )  POT 
u  =  (  I  )  A  T 
H  =11) 33 j M 
1  XI)  M  J  T 
0  <  =  <  I  )  U  3 
-XX ( I )SHIU3 
Jd= ( I ) s jMD 

i  a- ( ) x  p a s ; 

He = ( I ) 33 
0  = ( I ) SHAD 
se=(j)()A0j 
o  x  I ) r-  od 

V  XI)  H  3 
:0<X  i  )  ar 
P  X  1  )  DYM 
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FUA (  I  ) =SS 
FA (  I  ) =ST 
Z ( I ) =SU 
E ( I ) =SV 
AF (  I  )  =  SW 
C02 (  I  )  =SX 
02(1)  =S Y 
CO ( I ) =5Z 
LB (  I  )  =  3Z  A 
KN  T=  I 

C  KNT  GIVES  LOCATION  OF  LAST  INTERCHANGE 
C  SUCCEEDING  VALUES  IN  ORDER 

15  CONTINUE 
IF(KNT.EQ.l)  GO  TO  16 

C  IF  KNT  =  1  ALL  RATtOS  IN  ORDER 
LI M=KNT-1 
GO  TO  14 

16  CONTINUE 
C 

C  **  TO  RR I U  RESULTS  ** 

WR I T E ( 6  » 1 1 4 ) 

WR I TE ( 6  *  1 1 5  ) 

DO  17  1=1, NO 

WRITE(6»116)RUN( I),RATlu(l),RPM(l) , SPARK (  I ) , TOR ( I )  ,BHP(  I ) *  FhP ( I ) »  I 
1HP(  I  )  ,CF(  I ) » CBHP (  I ) ,MEFF(  I )  , LB { I ) , FU I  HP ( I ) , FU ( I ) ,CMEP( I ) ,ETA( I )  ,TE 
2M(  I  ) 

17  CONTINUE 

WR  I  TE ( 6 ,1 17 ) 

DO  18  1  =  1,  NO 

WRITE(6»118) RUN ( I ) , R A  T I g ( I ) , HYC ( I  ) >FA(  1  )  ,  FU A ( I )  ,Z(  I  )  »E(  I  )  ,  AF ( I )  , Cg 

1 M (  I  )  , MEN (  I  )  ,TS(  I  ) 

18  CONTINUE 

C  ***  FOR  PLOTTING  RESULTS*** 

IF(N0.LT.5)  GO  TO  24 
NN  =  NO 

C  TO  PLOT  CBHP,  FHP  AND  IHP 
RAT=RAT 10(1) 

L  K  =  0 
D  I  F  =  2  0  • 

T  =  0  •  0 

WR  I  TE ( 6  *  1 1 9 ) 

DO  19  1=1, NN 

CALL  PLOTER  ( NAT  I  0 (  I  )  ,CBHP(  I  )  ,FHP(  i  )  , i HP (  I  )  , T , D I F , RAT , LK ♦ ST^ I C , DD ) 

19  CONTINUE 

C  TO  PLOT  02*02  AND  CO 
RAT  =  RAT 10 (  1 ) 

L  K  =  0 
D  I  F  =  5  •  J 
T=  0  #  0 

WRITE(6,120) 

DO  20  1=1, NN 

CALL  PLOTER  (RAT  ILK  I  )  ,02 ( I ) , 02 ( I ) , CO ( 1 )  , T , D I F , RAT , LK , STO I C  ,  DD  ) 

20  CONTINJE 


)  :  i  .  » 

i  J  M  T  Z 


*1  » 


a-T. 


JVl  t  i  V 


C  =  (  I  )  A  a  1 
'!£  =  (!  )  A  3 
00  =  (  I  )  3 

ve= (i)3 

A  =  (  I  )  3  A 
X0* ( I ) S03 
Y  -  (  1  )  S  v 
Sc  =  (  I  )  03 
ASc  =  (  I  )  8J 
1  =TZ> 

j<_a' AH3H3  T *i  1  T  AJ  30  KOITA3QJ  23VI0  T/X 
1  J0  Hi  .  IUJA  J  uZ  I  Cl  3  11  J  2 

3U/ITM03  ei 
d  1  ,  T  00  (  I  #  L  1 .  T  Z  X  )  3  I 

/i  (  J T  T  A H  JJA  IsT/a  31 
i-i  /.**(•  i  j 
M  OT  00 
3UK.1T/iO0  di 


*  *  1  -iUOjM  T*  1  fls  OT  ** 

(  A I  I  *<3  )  JTlflW 

iexr*d)  jt la* 

0/  •  I  =  I  T I  OCT 

(  i  )  kiH<. »  l  :  j  ,1/.,  m  *  (  1  )  i*.  -i/i «  (  1  )  w  I  T  A>i  ♦  II  )  ht/i  (  d  I  *  d  )  -  T  13a 

"  1  1  >  •  1  O  •  l  *  (  l  )  u  j  *  (  1  )  -i  i  *  (  i  )  a  a3  *  (  1  )  3  3  ♦  (  i  )  Sr.  j 

(  I  )  MS 

30/  I  T  /l 03  VI 
(  V  1 1 ,  a  )  3  T  I  9\  a 
OK* 1=1  81  Ou 

<  i  )  J  *  '  .  )  3*  (  1  )  h  j'i  *  (  i  )  A  i»tl)3Vn#(I)^j,T  A>,  f  (  A  )  „0/i  (  8  1 1  ♦  c  )  3  T  1  X .. 

(  1  )  1  «  Cl)Z3f^  ♦  {  1  )  r 

JU/IT/03  81 
***2TJU.  3>i  0/ 1  T  TO.  M  3  ***  3 

AS  OT  00  (  e  •  T ,  ,-/  )  31 

0/  =  // 

3M1  (I/A  AH3  *  hh  3  TO  JH  ;T 
( 1)01 TAfl=TAfl 
C  =  >  J 
•  v,  S  =  3  I  <J 
U  •  w  =  T 
(  e  1 1 «  6  )  .3  T  I  fl  W 
HZ*  1  =  1  91  00 

, ,  A/<  *  ii(i»  i  ,  (  i  )  am  «  (  i  )  m  *  (.1  )Snd3t  (  i  )«i  TAJ  J  ZJ  1  u.3  JjAj 

3U/1T/03  91 

J3  UJ  A  So  *SC  3  To  JH  OT 

( non  a3  =  t  ah 

J 

i  .e  =  3  iu 

0.v«T 

(oSItc)3Tlfl'A 
/*»'  •  I*  1  wS  00 

/•*  lid*  i  *  (  i  )o3*  (  i  )  (  i  )  t  i  )ol  i  ,)  »j[WjH  j  _j  a  3 

31 /IT/03  s 
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C  TO  PLOT  TH.  E FF • »  TORQUE  AND  LOAD 
RAT=RA F I 0 ( 1 ) 

LK  =  0 
DI  F  =  2.  a 
T=0. 0 

WRI  TE( 3,121 ) 

DO  21  1=1 >NN 

CAi-L  FoOTE*  (  RAT  I  O  (  I  )  » ETA  (  I  )  ,ToR(  I  )  ,  Load  (  I  )  ,  T,DIF,*AT ,Lk,stoic,DD) 

21  CONTINUE 

C  TO  PLOT  MECH.  EFF * , C YL I NDER  TEMP./10.0  ,  AND  CBMEP 
rat  =  rat IOC  1 ) 

LK  =  0 
DI F  =  0.8 
T  =  2  0  •  0 

WRlTE(6,122) 

DO  22  I=1,NN 
C Y  T  EM (  I  )  =  TEM ( I ) / 1 0 • 0 

CALL  p  -UTER  (RATIO  (  I  )  ,MEFF  (  I  )  *  CYTEM  (  I  )  >CP'ER  (I),T,DIF,KAT»l-n,ST^IC 
1  » DD  ) 

22  CONTINUE 

C  TO  PLOT  FUEL  PER  CBHP»FUEL  PER  IHP  A*mDA/F  FRoM  ORSAT  ANALYSIS 
RAT  =  RAT 10 ( 1 ) 

LK  =  0 

D I  F= 10  1  • 

T=  2  0  •  0 

WRI TE( 3*123) 

DO  23  1=1, NN 
AIR (  I  )  =  10 • 0*F A ( I ) 

CALL  PLOTER  (RATIO!  I  )  *  FU  (  I  )  ,  FU  I  HP  (  1  )  ,  A 1 R  (  I  )  ,  T  ,  D  i  F  ,  RAT  ,  LK. ,  ST^  I  C  ,  Du  ) 

23  CONTINUE 

24  CONTINUE 

99  FORMAT! 13, A6 ) 

100  FORMAT ( F7. 3,F6.2 ,2  I  3, F4. 1 ,3A6 ) 

101  FORMAT ( 1HJ*F7.3»F6.2»2I3,3A6) 

102  FORMAT ( 1HJ ,68HRUN  SPARK  REVS  TIME  AIR  DP  CY  TEM  CU2  02  Co 
1  LOAD  COMMENTS) 

103  FORMAT ( 1HJ  ,79HRUN  SPARK  REVS  TIME  FUTEM  PRES  AIR  DP  CY  TEM  Co 

12  02  CO  LOAD  COMMENTS) 

104  FORMAT ( 1H1  ,40X,17H  CFR  TEST  RESULTS) 

105  FORMAT ( 1HL  ,  10X  ,9HTtST  NO*  ,I3,5X,6H  C * R .  =  ,  F6 • 2 , 5X , 9H  FULL  IS  ,*A6, 
14X,22H  BAROMETRIC  PRESSURE  =  ,  F7. 3 , 1  OX  ,  1 3 ) 

106  FORMAT (  1 2  ,  I  3 » F6 • 0 , F6 . 3 , F4 . 0 , F5 . 2 , F4 . 0 , F6 • 3 , F6 . 0 , F 5 . 1 , 2F4 . 1 , F5 . 1 , 3 A 
16, A4  ) 

111  FORMAT (  1HJ,  I3,I5,1X,F6.C*F6.3,F4.0,F6.3»F6.0,F6.1,2F5.1,F6.1 , 3A6 ) 

112  FORMAT ( 1H J  ,  I  3 , I5,1X,F6.0,F6.3,F5.0,F6.2,F4.0,F6.3,F6*0»F6.1,2F5.1, 
1F6. 1 »  3 A6 *  A4 ) 

113  FORMAT ( 1HJ»15X»14HFRICTI0N  LOAD= , F5 • 1 , 5H  LBS.) 

114  FORMAT( 1HK,40X,13H**  RESULTS  ** ) 

115  FORMAT ( 1HJ *25HRUN  F/A(MEAS  RPM  SPARK , 2X , 45HTORQUE  BHP  FH 

IP  IHP  CF  CBHP»2X»  8HMECH  EFF , 2X, 21HFU/HR  FU/IHP  FU/CoH 

2P,3X,23HIMEP  THER  EFF  CYL  TEMP) 

116  FORMAT ( 1HJ,I3,F8*4,F8*1,I4*2X,6F8*3,F7*2,3X,2F6»3,F8*3,3X,F  I  *  2 , F  8 • 
12,3X,F6*0) 
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0  =  >J 
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117  FORMAT ( 1HK  >  5 1H  RUN  F/A  MEAS  H/C  F/A  CAL  F/A  CAL  MOLES  ERRuK, 
HUH  A/F  RATIO, 8X, 8HC0MMENTS > 

118  FORMAT (1HJ>I3,5F8.4,2F8.2,5X,3A6»A4) 

119  FORMAT ( 1H1 > 10X  ,  38HCBHR , FHR , AND  I  HR  VERSUS  FUEL  Aik  RATIO) 

120  FORMAT ( 1H1 ,10X,37H  CU2 »  w2  AND  Co  VERSUS  FUEL  AIR  RATIO) 

121  FORMAT ( 1H1  ,1UX,53H  THERMAL  EFF.,  ToROUE  AND  LOAD  VERSUS  FUEL  AIR 
1RAT 10) 

122  FORMAT ( 1H1 »  10X  »  5  3H  MECH  EFF,  CYL  TEMP/10  AND  CBmER  VERSUS  FUEL  AIn 
1RAT 10) 

123  FORMAT ( 1H1  , 10X ,67HFUEL  PER  CBHR ,  FUEL  PER  I  HR  AND  FUEL  AIR  RATIo  F 
1R0M  ORSAT  ANALYSIS  ) 

STOP 

END 
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SIBFTC  PLOTER  DECK 

SUBROUTINE  PLOTER(X,  Y  ,R>S»T  ,DIF,RAT  , LK,sT^IC,DD  ) 
DIMENSION  CHAR (99) 

DATA  BLANK, PL0TY*PL0TR,PLuTS/1H  * 1H*, 1H+, 1HX/ 

IF ( LK.GT.  3  )  G  )  TO  1 

IF ( RAT  ,GT. STOIC )  GO  TO  1 

WRITE (6 *7)  STOIC 

1  CONTINJE 
LK=  1 

I F ( RAT.LE.STOTC )  GO  TO  5 
IF ( X.GT. STOIC'  GO  TO  5 
KK=( RA r~S  TO  I C ) *DD+ • 5 
DO  2  L= 1 »  K K 

2  WR I TE (  S  *  3 ) 

W  R  I  T  E  (  6  »  7 )  STOIC 
RA  T  =  ST  3  I C 

5  CONTINUE 

DO  6  JK  =  1 »  99 

6  CHAR ( JK ) =BLANK 
KK= ( RA  T-X ) *DD+ • 5 
I=Y*DIF~T 
J=R*DIF~T 
K=S*DIF-T 

IF ( I .LE.99 )  GO  TO  10 
1  =  99 

10  IF ( I • GT • 0 )  GO  TO  11 
1  =  1 

11  I F ( J  •  LE.99)  GO  TO  12 
J=  99 

12  I F ( J.Gf.O)  GO  TO  13 
J=1 

13  IF(K. LE.99)  GO  TO  14 
K=  9  9 

14  IF ( K.GT.O )  GO  TO  15 
K=  1 

15  CONTINUE 

I F ( J.EQ.I )  J=I+1 
IF(K.EG.J)  K= J+ 1 
IF(K.EQ.I)  K=  I  + 1 
CHAR (  I  )  =PLOTY 
CHAR ( J) =PLOTR 
CHAR ( K ) =P  LOT  S 
DO  8  L= 1 »  KK 

8  WR I TE ( 6  »  2 ) 

WR  I  TE  (  i>,4)  X,  Y,R,S,CHAR 
RAT  =  X 

3  FORMAT ( 1H J ) 

4  FORMAT ( 1H  » F 8 . 4 » 3 F7 . 2 » 1 X , 1 H . , 9 9A 1 ) 

7  FORMAT ( 1HJ »F8.4*23H  STOICHIUMtTRlC  MI XTURt *  10 ( 1H. * 9X)  ) 

9  RETURN 
END 
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APPENDIX 


SAMPLE  SET  OF  RESULTS 


C  F  R  f  t  S  1  RESULTS 


i/1 

■JJ 

h~ 

1  NO. 

5 

C  .  R  . 

Q  P  ? 

FUEL 

IS  PROPANE 

it- 

i E  S  U  L  i  S 

js-  «• 

R  U  ■  v 

F/ A ( ME AS 

R  PM 

SPARK 

1  Ok  (POE 

!  ’)  H  p 

FHP 

IMP  C 

F 

48 

0.05 74 

a  n  ^  *- 

9 

20.865 

2.  396 

0.977 

3.662  1 

.085 

4  9 

0.0501 

6 0  4*.  9 

21.829 

2.514 

0.979 

3.797  i 

.087 

80 

C  .0413 

594 . 4 

J 

20.5  14 

2.321 

0.962 

3.569  i 

.087 

D  1 

0-0320 

60  0. 1 

10 

15.3 4  2 

1.  7  93 

0.97  2 

2.964  1 

.08  8 

9  2 

0.0254 

596.6 

2  0 

8 » 0  6  5 

0.916 

0.966 

2.048  1 

.  08  8 

CBFiP 

4  £CH  E  F  F 

F  <j  /  H  R 

FU/IHP 

FU/CbHP 

IrtEP 

THEK  E  F  F 

SYL  TE 

2 . 6  8  5 

73.31 

1.670 

0.456 

0 . 6ZA 

128.71 

27.99 

8  0  3 

2.818 

74.21 

1.448 

0.381 

0.514 

133.18 

33.50 

856 

2 «  60  7 

73.04 

1  .  193 

0.334 

0.45  7 

127.41 

3  o  .  2  2 

8  30 

1.992 

67 .22 

0.938 

0.316 

0.471 

104.79 

40.36 

703 

1.082 

52.83 

0.7  53 

0.368 

0.696 

72.82 

34.71 

53  5 

